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ABSTRACT 


Ad  analytic  lolutlon  U  obtained  for  the  acouitic  prenure  statistics  in  a  closed  rectangular  shaped  cavity 
behind  a  simply  supported,  rectangular  plate  excited  by  boundary  layer  turbulence.  The  contribution  of  the 
cavity  acoustic  pressure  is  neglected  as  contributing  tc  the  plate  excitation,  leaving  only  the  turbulent  pres¬ 
sure  fluctuations  as  the  exciting  force.  The  mathematical  model  for  the  turbulent  pressure  statistics  is  based 
on  that  of  Corcot,  which  agrees  well  with  experiment.  A  byproduct  of  this  analysis  it  an  analytic  solution  for 
the  turbulent  flow  excited  plate  vibration  velocity  statistics.  The  plate  velocity  and  cavity  acoustic  pressure 
statistics  are  expressed  in  the  form  of  cross  power  spectral  densities  and  power  spectral  densities.  Dimension¬ 
less  forms  of  the  plate  velocity  spectral  density  and  cavity  acoustic  pressure  spectral  density  are  developed. 

The  dimensionless  plate  velocity  spectral  density  and  dimensionless  cavity  acoustic  pressure  spectral  den¬ 
sity  were  computed,  by  means  of  a  digital  computer,  for  selected  values  of  dimensionless  input  parameters. 
From  these  computed  dimensionless  spectra,  the  effects  of  major  parameters  on  the  plate  velocity  spectral 
density  and  the  cavity  acoustic  pressure  spectral  density  were  determined. 

A  'peak  spectrum,*  constructed  by  connecting  the  major  spectral  peaks  in  the  plate  velocity  or  cavity 
acoustic  pressure  .^ectra,  proved  to  be  a  useful  engineering  concept.  Knowledge  of  the  'peak  spectrum*  it 
equivalent  to  knowledge  of  the  maximum  plate  velocity  or  cavity  acoustic  pressure  spectral  levels  for  a  par¬ 
ticular  set  of  input  parameters.  Bated  on  the  computed  dimensionless  spectra,  mathematical  expressions  are 
derived  for  the  dimensionless  plate  velocity  'peak  spectral  density*  and  the  cavity  acoustic  prenure  'peak 
spectral  density*  over  a  limited  range  of  dimensionless  frequency.  The  computed  simply  supported  plate  veloc¬ 
ity  'peak  spectrum*  compares  well  with  the  plate  velocity  'peak  spectrum*  constructed  from  experimental 
measurements  on  a  fixed  edge  plate  above  the  first  plate  natural  frequency.  No  experimental  data  exists  for 
the  cavity  acoustic  pressure. 

Comparison  of  the  computed  dimensionless  cavity  acoustic  pressure  spectral  density  at  the  plate  and  the 
dimensionless  turbulent  pressure  spectral  density  allowed  formulation  of  criteria  under  which  the  cavity  acous¬ 
tic  pressure  wasoegliglble  compared  to  the  turbulent  pressure.  As  this  analysis  assumed  the  cavity  acoustic 
pressure  to  be  negligible  compared  to  the  turbulent  pressure,  the  aforementioned  criteria  are,  in  effect,  limits 
of  applicability  of  this  analysis. 
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I 

INTRODUCTION 

Flow  induced  noise  is  largely  responsible  for  limiting  submarine  and 
surface  ship  sonar  performance,  producing  objectionable  noise  levels  Inside 
high  speed  aircraft,  and  producing  relatively  high  vibration  levels  in  the  pro¬ 
pulsion  systems  of  modern  spacecraft  which  increase  the  probability  of 
failure  of  system  components.  Because  of  this  wide  variety  of  problems, 
theoretical  and  experimental  research  in  flow  noise  and  flow  induced  noiBe  has 
increased  in  the  past  decade.  The  research  in  flow  noise  has  been  aimed  at 
defining  a  statistical  model  for  the  turbulent  boundary  layer  pressure  and/or 
velocity  which  provides  the  excitation  to  the  mechanical  system.  The  research 
in  flow  induced  noise  has  been  primarily  aimed  at  theoretical  solutions  for  the 
response  characteristics  of  various  systems  using  mathematically  tractable 
approximations  to  the  turbulent  boundary  layer  excitation.  Certain  experimen¬ 
tal  studies  of  flow  induced  noise  in  simple  systems  have  also  been  performed. 

To  date,  although  there  is  general  agreement  in  the  measurements  of 
boundary  layer  pressure  s’atistics,  there  is  not  complete  agreement  as  to  the 
mathematical  model  of  these  statistics.  Also,  although  a  wide  variety  of  ele¬ 
mentary  flow  induced  noise  problems  have  been  studied,  few  have  used  any  of 
the  existing,  experimentally  based,  mathematical  models  for  the  boundary 
layer  excitation. 
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2 


It  is  the  purpose  of  this  study  to  add  to  the  understanding  of  flow  induced 
noise  by  investigating  the  effects  of  major  parameters  on  the  sound  field  pro¬ 
duced  in  a  closed  space  behind  a  simply  supported  plate  excited  by  boundary 
layer  turbulence  utilizing  an  existing,  experimentally  based  mathematical 
model  of  the  turbulent  wall  pressure. 


II 


OBJECTIVES 

The  objective  of  this  study  is  to  provide  the  submarine  sonar  systems 
designer  with  information  concerning:  (1)  the  acoustic  environment  of  sonar 
transducers:  and  (2)  the  major  parameters  which  may  affect  this  environment. 

The  above  objective  has  been  attacked  by  means  of  an  analytical  study  of 
the  acoustic  field  in  a  closed  space  behind  a  simply  supported,  rectangular, 
flat  plate  which  is  excited  by  turbulent  boundary  layer  pressure  fluctuations. 
The  closed  space  is  bounded  by  five  rigid  walls  and  the  flexible  plate.  A 
sketch  of  the  model  used  in  this  analysis  is  shown  in  Figure  1.  This  model 
provides  a  fair  representation  of  the  acoustic  environment  of  sonar  trans¬ 
ducers  in  submarines. 
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TECHNICAL  APPROACH 


3. 1  Technical  Background 

The  purpose  of  this  section  is  twoiold.  First,  it  summarizes  work  in  the 
field  of  flow  and  flow  induced  noise.  Secondly,  it  provides  justification  for 
certain  assumptions  made  in  the  analysis  to  follow. 

3. 1. 1  Theoretical  Studies  in  Flow  Induced  Noise 

Early  theoretical  studies  in  the  field  of  flow  induced  noise  were  directed 
toward  a  prediction  of  noise  in  aircraft  fuselages  [l,  2].  The  mathematical 
models  of  the  turbulent  boundary  layer  in  these  studies  were  not  based  on 
experimental  evidence.  During  this  same  period,  Lyon  [3]  studied  the  re¬ 
sponse  of  strings  to  random  excitation,  Eringen  [4]  derived  expressions  for 
the  response  of  beams  and  plates  to  random  pressure  fields,  and  Kraichnan  [5] 
studied  the  free  radiation  of  sound  from  turbulent  excitation  of  a  series  of  thin, 
stiff  flat  plates.  Again,  the  excitation,  although  characteristic  of  boundary 
layer  turbulence  in  certain  respects,  was  not  compatible  with  experimental 
data. 

Dyer  [fi]  was  one  of  the  first  to  study  the  coupled  plate  vibration-acoustic 
radiation  i  roi.Jcm.  He  assumed  the  boundary  layer  pressure  correlation  func¬ 
tion  to  be  In  .'  product  of  a  convected  spacial  delta  function,  a  fixed  spacial 
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delta  function,  an  amplitude,  and  a  decaying  function  of  time.  Although  the 
delta  functions  were  relatively  poor  approximations  to  the  actual  longitudinal 
and  lateral  pressure  correlations,  the  convection  and  time  decay  were  phe¬ 
nomena  which  agreed  with  experiments.  The  above  excitation  was  used  in  a 
normal  mode  approach  to  the  prediction  of  the  acoustic  field  in  a  closed  space 
behind  a  simply  supported,  flow  excited,  flat  plate.  The  walls  enclosing  the 
space  were  pressure  release  surfaces.  Although  Dyer's  input  waa*not  precise, 
it  provided  insight  into  the  behavior  of  the  model  considered  and  encouraged 
further  work  utilizing  normal  mode  theory. 

A  short  time  layer,  Dyer  [7]  used  the  same  model  of  the  turbulent  pres¬ 
sure  correlation  to  calculate  the  displacement  correlation  function  of  a  turbu¬ 
lence  excited  flat  plate.  About  this  same  time,  Strasberg  [8]  used  a  mathe¬ 
matical  model  of  the  turbulent  pressure  correlation  based  on  the  data  of 
Harrison  [9]  to  predict  the  displacement  spectral  density  of  plates  and  mem- 
brances.  This  was  one  of  the  first  cases  where  experimentally  based  inputs 
were  used. 

Powell  [10]  investigated  the  fatigue  of  structures  excited  by  random  pres¬ 
sure  fields  and  indicated  that  the  response  cross  spectral  density  is  maximized 
when  the  incident  pressure  correlation  matches  the  modal  wavelength.  Al¬ 
though  this  is  a  possible  condition  in  spacecraft  because  of  their  high  speeds, 
it  is  not  likely  to  occur  in  submarines. 

Maidanik  and  Lyon  ["]  studied  the  response  of  strings  to  moving  noise 
fields  using  the  Dyer  delta  function  model  for  the  pressure  correlation. 


From  1960  to  the  present,  an  extensive  research  program  in  flow  noise 
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and  flow  induced  noise  has  been  conducted  at  the  University  of  Southampton. 
This  program  yielded  a  recent  paper  by  Mercer  [12]  in  the  response  of  multi- 
supported  beams  to  a  random  pressure  field. 

Tack  and  Lambert  [13]  derived  general  expressions  for  the  response  of 
plates  and  bars  to  boundary  layer  turbulence. 

A  recent  paper  by  White  [14]  used  an  experimentally  based  expression 
for  the  turbulent  pressure  cross  spectral  density  to  predict  the  response  and 
consequent  sound  radiation  from  a  rectangular  flat  plate.  The  response  of  the 
plate  and  the  acoustic  radiation  are  averaged  over  frequency  bands  bo  that  the 
details  of  the  response  are  bypassed,  and  only  the  effect  of  such  gross  param¬ 
eters  as  panel  modal  density  and  boundary  layer  characteristics  are  investi¬ 
gated.  White  obtains  very  good  agreement  between  theory  and  experiment 
within  these  assumptions. 

Bull  et  al.  [15],  as  a  part  of  the  University  of  Southampton  effort,  calcu¬ 
lated  the  displacement  spectral  density  of  flat  plates  due  to  turbulent  excitation. 
They  used  a  normal  mode  approach  to  the  plate  problem  and  an  experimentally 
based  expression  for  the  turbulent  pressure  correlation  as  input.  Collier  [l6], 
using  an  approximation  to  the  turbulent  pressure  correlation  function,  calcu¬ 
lated  the  plate  acceleration  correlation  functions  and  acoustic  pressure  field 
radiated  from  the  plate  into  an  infinite  fluid  space  at  rest. 

Pretlove  [17]  presents  the  theory,  from  a  normal  mode  approach,  for 
calculating  the  displacement  spectral  density  of  a  simply  supported  panel, 
backed  by  a  rectangular  closed  cavity.  In  the  forcing  function  for  the  plate,  he 
makes  provision  for  both  the  turbulent  pressure  and  the  resultant  acoustic 
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pressure  in  the  cavity.  He  does  not,  however,  solve  this  problem  for  any 
particular  model  of  the  turbulent  pressure.  From  his  theory,  however,  he 
states  that  the  effect  of  the  backing  cavity  is  most  severe  in  the  cases  of  thin 
panels  covering  shallow  cavities. 

3. 1.  2  Experimental  Studies  in  Flow  Induced  Noise 

Few  experimental  studies  of  flow  Induced  noise  appear  in  the  literature. 

El  Baroudi,  Ludwig,  and  Rlbner  [is]  experimentally  investigated  the  displace¬ 
ment  correlation  properties  of  a  flow  excited  plate  and  the  resultant  total  sound 
power  radiated  into  a  reverberant  room.  In  a  related  effort,  el  Baroudi  [19] 
measured  the  displacement  correlation  properties  and  displacement  spectral 
density  of  thin  flat  plates  excited  by  turbulent  boundary  layers. 

Bull  et  al.  [15]  measured  the  displacement  spectra  of  thin  flat  plates  due 
to  turbulent  excitation.  The  plates  were  effectively  mounted  in  a  fixed  manner, 
and  the  experimental  results  were  compared  to  the  theoretically  predicted 
spectra  for  the  simply  supported  case  using  an  experimentally  based  model  for 
pressure  correlation  as  input.  As  might  be  expected,  agreement  between 
theoretical  and  experimental  results  was  not  good. 

Maestrello  [20]  measured  the  sound  power  spectra  radiated  into  a 
reverberation  chamber  from  turbulence-excited  plates.  In  another  paper  [2l] 
he  measured  the  plate  displacement  correlation  function  for  turbulence  excited 
plates. 

3. 1. 3  Theoretical  Studies  of  Flow  Noise 

Flow  noise,  as  used  in  this  study,  refers  to  the  wall  pressure  fluctuations 
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produced  by  a  turbulent  boundary  layer.  The  mathematical  theory  of  pressure 
fluctuations  in  homogeneous,  isotropic  turbulence  is  well  developed  [22,  23, 

24,  25,  26],  However,  because  of  its  extreme  complexity,  the  problem  of 
pressure  fluctuations  produced  by  a  turbulent  boundary  layer  has  not  been 
treated  extensively. 

Kraichnan  [27,  28],  Lllley  and  Hodgson  [29] ,  and  Sternberg  [30]  have 
made  theoretical  studies  aimed  at  computing  the  mean  of  the  turbulent  wall 
pressure,  but  the  space-time  correlation  properties  (or  their  Fourier  trans¬ 
form)  are  required  in  order  to  treat  flow  Induced  noise  problems. 

Two  attempts  have  been  made  to  predict  the  correlation  properties  of 
turbulent  wall  pressure  fluctuations.  Gardner  [3l]  ,  starting  with  the  Navier- 
Stokes  equations,  attempted  to  predict  the  wall  pressure  correlation  by 
assuming  forms  for  the  various  velocity  correlations  occurring  in  his  expres¬ 
sions.  His  results  were  in  sharp  disagreement  with  experiment.  White  [32], 
claiming  that  Gardner's  approach  was  sound  and  that  mathematical  errors 
were  the  cause  of  the  discrepancy  between  theory  and  experiment,  attacked  the 
same  problem  through  Gardner's  approach.  White  calculated  the  turbulent 
wall  pressure  cross  correlation  function,  lateral  and  longitudinal  cross  spec¬ 
tral  densities,  and  convection  velocity.  White's  predictions  agree  reasonably 
well  with  experiment  except  in  the  case  of  the  wall  spectral  density  where 
there  is  a  large  discrepancy  between  predicted  and  experimental  results. 

White  shows  that  the  pressure  cross  spectral  density  is  approximately  equal  to 
the  product  of  the  lateral  and  longitudinal  cross  spectral  densities,  which 
bears  out  a  previous  empirical  prediction  by  Corcos  [33], 
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3. 1. 4  Experimental  Studies  in  Flow  Noise 

The  first  reliable  measurements  of  turbulent  boundary  layer  pressure 
statistics  were  made  by  Willmarth  [34]  in  1956.  Harrison  [9]  in  1958  pub¬ 
lished  the  first  measurements  of  the  wall  pressure  cross  spectral  density. 
Since  these  measurements,  many  experimenters  have  measured  the  turbulent 
wall  pressure  cross  spectral  density  and  correlation  functions.  Extensive 
bibliographies  of  this  work  are  presented  by  White  [35]  and  Bull  et  al.  [15], 
and  only  selected  references  will  be  presented  in  this  study. 

Bakewell  et  al.  [36]published  extensive  data  of  mean  square  pressure, 
pressure  spectral  densities,  and  pressure  correlation  functions  from  experi¬ 
ments  in  a  3-1/2-inch-diameter  turbulent  air  flow  facility. 

Using  both  theoretical  arguments  and  data  obtained  at  the  Ordnance 
Research  Laboratory  at  Pennsylvania  State  University,  Skudrzyk  and  Haddle 
[37]  derived  the  following  semi-empirical  expression  for  the  turbulent  wall 
pressure  spectrum  as  measured  with  very  small  transducers: 


=  0.75  x  io*vP(Ju»«r 
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where  a  is  a  constant  which  takes  on  different  values  for  different  fluids. 
For  water, 


and  for  air, 


a  -  1.0 
a  -  }.0 


(3.  2) 


Although  this  expression  is  not  in  exact  agreement  with  experimental  data 


in  shape  because  of  the  discontinuous  nature  of  (3. 1),  it  shows  good  agreement 
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with  experimental  data  in  both  amplitude  and  shape  at  frequencies  above  the 
discontinuity  and  a  reasonable  approximation  to  the  amplitude  at  frequencies 
below  the  discontinuity.  The  agreement  between  (3. 1)  and  experiment  1  data 
will  be  discussed  further  in  Section  4. 1. 

Wlllmarth  and  Wooldridge  [38]  reported  further  measurements  in  1962, 
designed  to  provide  more  detailed  statistical  properties  of  the  wall  pressure 
fields.  About  this  same  time,  Serafini  [39]  published  new  data  of  wall  pres¬ 
sure  measurements. 

Corcos  [33,  40],  using  the  data  of  Wlllmarth,  Bakewell,  and  Serafini, 
stated  that  the  turbulent  wall  pressure  cross  spectral  density  function  could  be 
expressed  as 

lV.  (3.3) 

Plots  of  A(<u£/Uc)  and  b(<j»;/Uc)  are  also  included  in  Corcos' work  and 
are  reproduced  in  Figures  2  and  3.  Corcos  went  on  to  predict  the  error  re¬ 
sulting  from  the  measurement  of  the  turbulent  wall  pressure  field  with  finite- 
size  transducers  and  presented  a  means  of  correcting  measurements  for  this 
error. 

Researchers  at  the  University  of  Southampton  have  made  extensive  meas¬ 
urements  of  the  turbulent  wall  pressure  statistics.  The  resulting  reports  are 
summarized  in  the  1963  report  by  Bull  et  al.[l5],  Wlllmarth  and  Roos  [4l] , 
incorporating  Bull's  data,  reattacked  the  problem  of  resolution  by  finite-size 
transducers  on  the  basis  that  Corcos'  similarity  form  for  the  cross  spectral 
density,  equation  (3.3),  although  accurate  over  a  wide  frequency  range,  failed 
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to  provide  accurate  hydrophone  size  corrections  at  high  frequencies.  Will- 
marth  indicates  that  the  reason  for  this  limitation  is  that  (3.  3)  is  inaccurate  at 
small  spacial  separations.  However,  as  Willmarth  points  out,  wall  pressure 
correlation  measurements  at  small  spacial  separations  have  never  been  made 
because  of  experimental  difficulties.  He  therefore  does  not  propose  an  alter¬ 
nate  form  to  (3. 3). 

Some  of  the  data  mentioned  above  have  been  obtained  at  the  walls  of  circu¬ 
lar  pipes  and  others  at  the  surfaces  of  flat  plates.  The  lack  of  any  significant 
difference  between  the  pressure  spectra  obtained  in  both  cases  is  indicated  in 
a  recent  work  by  Schloemer  [42],  who  measured  turbulent  wall  pressure  data 
in  the  presence  of  favorable  and  adverse  pressure  gradients. 


IV 


THEORETICAL  DEVELOPMENT 

The  model  treated  in  this  analysis  was  a  simply  supported,  rectangular, 
flat  plate  mounted  in  an  infinite  rigid  baffle  and  backed  by  a  rectangular  cavity 
with  rigid  walls.  Fluid  flows  over  the  top  of  the  rigid  baffle  and  plate,  and  the 
cavity  is  filled  with  an  acoustic  fluid.  A  graphic  representation  of  this  model 
is  presented  in  Figure  1. 

In  order  to  determine  the  acoustic  pressure  field  within  the  cavity,  it 
would  ordinarily  be  necessary  to  determine  the  vibration  of  the  plate  due  to  the 
combined  turbulent  and  cavity  acoustic  pressure  excifations.  The  cavity  acous¬ 
tic  phase  velocity  at  the  plate  would  then  be  equated  to  the  plate  velocity,  re¬ 
sulting  in  an  integral  equation  for  the  acoustic  velocity  potential.  The  various 
constants  resulting  from  the  solution  of  this  integral  equation  would  be  deter¬ 
mined  by  the  remaining  boundary  conditions,  which  require  that  the  acoustic 
phase  velocity  be  zero  on  the  rigid  cavity  walls.  Pretlove  [17]  used  an 
approach  similar  to  the  above  in  computing  the  plate  displacement  spectral 
density  for  a  model  identical  to  that  described  above.  However,  Pretlove 
discovered  that  the  cavity  acoustic  pressure  had  little  effect  on  the  plate  dis¬ 
placement  spectral  density  except  in  cases  of  thin  plates  covering  shallow 
cavities.  For  submarine  applications,  the  plates  covering  the  cavities  are 
not  thin  and,  therefore,  in  accordance  with  Pretlove' s  results,  it  is  assumed 
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in  this  analysis  that  the  cavity  acoustic  pressure  provides  negligible  excita¬ 
tion  to  the  plate;  that  is,  the  only  forces  exciting  plate  vibrations  are  those 
associated  with  the  turbulent  boundary  layer  pressure  fluctuations. 

The  above  assumption  simplifies  the  problem  considerably.  To  compute 
the  cavity  acoustic  pressure  field,  it  is  now  merely  necessary  to  compute  the 
plate  vibration  velocity  due  to  the  turbulent  boundary  layer  excitation  and  to 
equate  this  result  to  the  cavity  acoustic  phase  velocity  at  the  plate.  The  re¬ 
sulting  acoustic  velocity  potential  must  also  satisfy  the  zero  velocity  condition 
at  the  rigid  walls.  Hence,  it  is  only  necessary  to  determine  a  model  for  the 
turbulent  boundary  layer  pressure  statistics  in  order  to  solve  the  above 
problem. 

Since  the  turbulent  boundary  layer  pressure  is  a  random  phenomenon, 
the  cavity  acoustic  pressure  resulting  from  the  boundary  layer  excitation  of 
the  plate  will  also  be  a  random  phenomenon.  Random  processes  are  usually 
treated  in  terms  of  their  space-time  correlations  or  in  terms  of  their  sped  ral 
properties  in  the  frequency  domain.  Collier's  recent  work  [l6]  on  the  vibra¬ 
tion  and  acoustic  radiation  of  turbulence  excited  plates  indicated  certain  com¬ 
putational  difficulties  arising  from  treatment  of  the  statistical  properties  in 
the  time  domain.  Therefore,  in  this  analysis,  random  variables  will  be 
described  in  terms  of  their  spectral  properties  in  the  frequency  domain. 

4. 1  Mathematical  Model  of  the  Turbulent  Wall  Pressure  Cross  Spectral  Density 

The  mathematical  model  of  the  turbulent  wall  pressure  cross  spectral 
density  selected  for  use  in  this  analysis  was  that  of  Corcos  [33]  ,  which  is 
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given  by  (3.3)  and  is  repeated  here  for  reference. 

l,<>.  (3.3) 

Corcos'  model  was  selected  for  two  reasons.  First,  its  mathematical 
form  is  such  that  space  variables  are  separated,  thereby  easing  computational 
difficulties.  Secondly,  Corcos'  model  was  used  to  theoretically  predict  trans¬ 
ducer  size  corrections  which  agree  well  with  measurements  over  a  wide  range 
of  frequencies.  Thus,  subject  to  the  limitation  at  high  frequencies  discussed 
in  Section  3. 1.4,  Corcos'  model  of  the  turbulent  wall  pressure  cross  spectral 
density  provided  a  realistic  and  mathematically  attractive  model  for  this 
analysis. 

Before  proceeding  further,  it  is  necessary  to  discuss  the  limitations  in 
applicability  of  equation  (3. 3).  Since  £  and  rj  are  relative  coordinates, 
equation  (3.3)  tacitly  assumes  homogeneous  stationary  turbulence.  Therefore, 
Corcos'  model  of  the  cross  spectral  density  is  strictly  applicable  only  to  the 
case  of  a  flow  having  zero  pressure  gradient  and  constant  boundary  layer 
thickness.  However,  equation  (3.  3)  can  be  used  with  good  accuracy  for  turbu¬ 
lent  flows  in  which  slow  boundary  layer  growth  and  small  pressure  gradients 
occur.  In  this  analysis,  the  flow  is  assumed  to  have  the  following  charac¬ 
teristics: 

a.  constant  boundary  layer  thickness  over  the  plate, 

b.  zero  pressure  gradient. 

Corcos  [33]  also  shows  that  the  convection  velocity  is  a  function  of  the  Strouhal 
Number  ,„8*  uo .  This  relationship  is  shown  in  Figure  4.  Since  the  variation 
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of  the  convection  velocity  is  not  large  over  a  wide  range  of  Strouhal  Number, 
it  is  further  assumed  in  this  analysis  that  the  convection  velocity  is  a  constant 
given  by 

U  -0.65U  (4.1) 

c  o 


Figures  2  and  3  present  curve  fits  to  experimental  data  for  the  functions 
A(kj£/Uc)  and  B(wj//Uc)  contained  in  (3. 3).  From  these  figures,  the  follow¬ 
ing  expressions  have  been  selected  to  represent  the  A  and  B  functions  based 
on  a  balance  of  curve  fit  and  mathematical  simplicity: 


and 
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(4.3) 


Note  that  equation  (4.  2)  fits  the  measured  value  of  a(w£/Uc)  very  well, 
whereas  equation  (4.3)  sacrifices  some  accuracy  in  fitting  the  experimental 
data  in  favor  of  mathematical  simplicity.  The  expression 
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provides  a  much  better  fit  to  the  experimental  data  than  does  equation  (4.3), 
but  its  use  in  the  analysis  to  follow  greatly  increased  the  mathematical  com¬ 
plexity.  Hence,  equation  (4.3)  was  selected  to  represent  B(o>»;/ut). 

It  remains  to  select  an  expression  for  <t>(cu)  in  equation  (3.  3)  in  order  to 
completely  specify  the  mathematical  form  of  the  turbulent  wall  pressure  cross 
spectral  density.  Lilley  and  Hodgson  [29] ,  and  Skudrzyk  and  Haddle  [37]  have 
proposed  expressions  designed  to  describe  the  turbulent  wall  pressure  spectral 


density.  The  agreement  of  the  Skudrzyk  and  Haddle  expression  with  experi- 
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mental  results  is  much  better  than  the  LUlcy  and  Hodgson  expression.  Fur¬ 
ther,  the  mathematical  form  of  the  Skudrzyk  and  Haddle  model  is  quite  simple. 
Hence,  the  Skudrzyk  and  Haddle  expression  for  the  turbulent  wail  pressure 
spectral  density  (equations  (3. 1)  and  (3.  2))  was  used  in  this  model  and  is  re¬ 
peated  here  for  reference: 

*  0.75  v  l(T'Vp(Ju’/>* 
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where 

a  =  1 .0  for  water 
a  =  3-0  for  air  . 

Figure  5  presents  a  comparison  of  the  experimental  data  of  Bakewell  [36]  , 
Schloemer  [42],  and  Bull  [15]  with  the  Skudrzyk  and  Haddle  model.  The  data 
are  lor  air  rather  than  water  because  of  the  greater  availability  of  air  data. 
The  agreement  between  the  experimental  data  and  the  Skudrzyk  and  Haddle 
expression  is  poor  at  low  non-dimensional  frequencies  but  improves  consid¬ 
erably  above  the  cutoff  frequency.  Although  the  data  of  Schloemer  and  Bull 
were  obtained  on  flat  plates  and  those  of  Bakewell  in  a  pipe  flow  facility,  the 
Skudrzyk  and  Haddle  non-dimensional  form  brings  these  data  into  excellent 
agreement.  It  may  therefore  be  assumed  that  the  parametric  form  of  equa¬ 
tion  (3. 1)  is  valid,  but  the  values  of  the  constants  and  cutoff  frequency  could 
be  considerably  improved.  For  this  analysis,  however,  equation  (3. 1)  was 
used  in  the  original  form  to  describe  the  turbulent  wall  pressure  spectral 
density. 

Combining  equations  (3.  3),  (4.  2),  (4.  3),  and  (3. 1),  the  turbulent  wall 


oi  <  1.256  — 

5* 


(3.1) 


U)  >  1.256  —  . 

FT 
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pressure  cross  spectral  density  may  be  described  by 


S|)pl  -  O.’S  >  I0‘  '  n3  f.j  ff  1'cXe'0-7la'T'/l,cl]ei,a'^  Uc> 


u 

oi  <  1.256  —  . 

"  5* 


2 

>  ({,n,oj)-  1.5  x  lo-’a2^.  [c-o.m|-f/utJe-o.7|a,vuc|]e-i^f/et) 
pp 


oi  >  1.256  — 

5* 


(4.4) 


Equation  (4.4)  is  the  mathematical  model  of  the  turbulent  wall  pressure  cross 
spectral  density  used  in  the  analysis  to  follow.  For  convenience  of  reference, 
the  assumptions  pertinent  to  equation  (4.  4)  are  summarized  as  follows: 

a.  constant  boundary  layer  thickness  over  the  plate, 

b.  small  pressure  gradients. 

It  is  further  assumed  that  the  ratio  of  the  convection  velocity  to  the  free 
stream  velocity  is  constant  as  presented  in  equation  (4. 1). 


4. 2  Development  of  the  Plate  Velocity  Cross  Spectral  Density 
4.  2. 1  Plate  Velocity  Response  to  a  Deterministic  Pressure 

The  differential  equation  governing  the  displacement  of  the  plate  due  to 
the  turbulent  boundary  layer  pressure  excitation  on  the  plate  surface  is 


nv 


4  dw 

w  4-  r  — 


dt 


(4.5) 


where 

D  is  the  flexuril  rigidity, 

r  is  the  effective  damping  coefficient  per  unit  area,  and 
p  is  the  effective  mass  per  unit  area. 

The  terms  "effective  mass"  and  "effective  damping"  as  used  above  denotes 

! 

; 


IB 


the  mass  and  damping  due  to  the  combined  effect  of  the  plate  and  water.  In 
this  treatment,  for  lack  of  detailed  information  concerning  the  effect  of  the 
water  on  mass  and  damping,  the  effective  mass  and  the  effective  damping 
coefficient  are  assumed  to  be  constants.  It  is  further  tacitly  assumed  in 
equation  (4.  5)  that  the  turbulent  pressure  field  is  not  affected  by  the  plate 
motion.  The  effects  of  cavity  acoustic  pressure  on  the  plate  and  neglected, 
and  the  forcing  function  is  only  the  turbulent  pressure. 

The  solution  to  equation  (4.  5)  for  any  arbitrary  deterministic  pressure 
field  can  be  determined  by  a  supeiposition  of  the  normal  modes  of  vibration 
of  the  corresponding  free-undamped  plate,  governed  by 

r>v  'w  ♦  [i  — -  -  o .  (4.6) 

Ot2 

The  solutions  to  equation  (4.  6)  satisfying  the  simply  supported  edge  con¬ 
ditions  shown  in  Figure  1  are  given  by 

w(x.y  .0  =  amn(i.y)  sin  wmnt,  (4.7) 


where  the  mode  shapes  and  corresponding  natural  frequencies  are  given  by 


and 


(4.8) 


(4.9) 


The  normal  modes,  «mn ,  form  a  complete  set  of  orthonormal  functions: 


that  is, 


a  (x  ,y) a  (x  ,v) dx  dy 
m  n  7  i|f  '  7 


=  8„s»  .• 


(4. 10) 


■  t  _  _ 
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It  is  now  assumed  that  the  solution  to  equation  (4.  5),  for  any  determinis¬ 
tic  pressure,  can  be  written  as  a  sum  of  the  normal  modes,  each  multiplied  by 
a  function  of  time:  that  is, 

■x. 

•v 

w(x  ,y  ,t)  -  u  (*,  y)T  It),  (4.  11; 

7  d—j  m  n  7  m  n 

•n  1 
n  =1 

where  T  is  to  be  determined. 

inn 

Substituting  (4. 11)  into  (4.  5)  and  utilizing  (4. 10),  one  finds  that  Tmn(t) 
must  satisfy 

dJT  dT  i  /b/i 

-£T  V—  +  p,(*.y.o «„,(*. rldidy  .  (4.12) 

At  this  point,  it  is  convenient  to  solve  equation  (4.  5)  by  meanB  of  equa¬ 
tions  (4. 11)  and  (4. 12)  for  two  special  cases.  For  the  first  case,  let 

P,(*.y.O  -*(*  -  *0 Sty  -y')*1"',  (4. 13) 


that  is,  a  concentrated  load  applied  at  (*',y')  varying  sinusoidally  in  time. 
Assume  the  solution  can  be  written  as 


(4. 14) 


Equation  (4. 14)  is  the  defining  statement  for  H(x,x\y,y',aj),  which  is  termed 
the  complex  frequency  response.  Solution  of  (4. 12)  with  p,(x,y ,  o  defined  by 
(4. 13)  results  in  the  following  solution  via  (4. 11): 


w(*,y,c)  - 


no 

E 

m  3  1 


XT7^7!£ 

L  M  . 


I«l 


(4.  15) 


Hence,  comparison  of  (4. 15)  and  (4. 14)  yields  the  following  solution  for  the 


i 


i 


i 
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complex  frequency  response: 


u,  /  ,  .  l  ^  a«#(* ly)0*#!*  >y ) 

H(«,«  ,y  ,y  ,w)  --  2_, 


*1  ma  ^ 


(4.16) 


As  the  second  special  case,  consider  the  impulsive  loading  at  (*',y‘) 
occurring  at  time  t  ': 


p,(*.y  ,0  -  5(*  -  i')5(y  -  y')  5(t  -  r#) 


(4. 17) 


and,  define: 


0  ■  t  - 1' 

/w  (*  ,  y  ,  t) 


h  (»  .  «',y  ,  y',  0) 


0 


0  >  0 


0  <  0  . 


(4.18) 


The  solution  of  equation  (4. 12),  when  pt  is  as  given  by  (4. 17),  is 


(4. 19) 


The  constants  C,  and  c2  may  be  evaluated  by  applying  initial  conditions 
appropriate  to  an  impulsive  loading,  with  the  following  result: 


T..w  -  i=z=  ...(«>•) ...  y-;.  -  (£j  e  . 

w 


Hence,  from  (4. 11)  and  (4.  18) 


-  f-fej 


(4.  20) 


(4.21) 


e  >o 


By  use  of  (4.  21)  and  the  principal  of  superposition  for  linear  systems,  the 
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plate  response  for  any  deterministic  pressure  excitation  may  be  constructed. 
This  is  accomplished  by  considering  the  excitation  pressure  field  to  be  the 
superposition  (or  summation)  of  an  infinite  number  of  impulses  in  time  and 
space.  Hence,  for  any  deterministic  pressure  field,  the  response  may  be 
written: 

w(*,y,c)»j^  p((*',y',t')h  (*,«',  y  ,y',0)dx'dy'dt\  (4.22) 

However,  as  0  =  t  -  t\  (4.  22)  may  be  rewritten  as 

—-nr  p^i'.y'.t-Olhii.i'.y.y',  0)dx'iy'iO  (4.  23) 


In  the  acoustic  problem,  the  velocity  of  the  plate  (rather  than  the  dis¬ 
placement)  will  be  required  for  the  boundary  value.  Hence,  defining 


<(*,*',y  ,y',  0)  - 


3h(«,i',y.y'.  0) 

dt 


(4.  24) 


as  the  velocity  response  of  the  plate  to  an  impulsive  loading,  from  (4. 24)  and 
(4. 21)  the  velocity  response  is  found  to  be 


d>  o . 

(4.  25) 


Defining  the  velocity  field  of  the  plate  as 


<?w(«,y  ,t) 


<£(*,y  .0 


dt 


(4.  26) 
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one  may  express  the  velocity  response  of  the  plate  to  a  deterministic  pressure 
excitation  via  (4.  25)  and  (4.  26)  as 

<£(«.y  .0  -  Ctf  P,(*'ty - 9)  £(« ,«',y  ,y\  9) dx' dy'  d  9 .  (4.  27) 


4.  2.  2  Plate  Velocity  Cross  Correlation 

The  plate  velocity  cross  correlation  is  defined  as 

Q^^(*r«2.y  ,.yj.‘,.,2)  “  E  l0(«|.y,.‘,>  ^(ij.y,.*,))  (4. 28) 

and  the  turbulent  wall  pressure  cross  correlation  as 

Qpp  (* i’*2'V i i •* 2^ "  EU>t<*,.y  |.*,>pt(*  jiy 


where  E  denotes  the  ensemble  average. 
From  (4.  27)  and  (4.  28), 


<W,i-,2,yi’yr‘i,l2>  ”  E 


p«<*>V^-*i>p«(x>2-‘2-fl2) 


ojojojojojo 


((*  i  >*',  <y  i  .y', .  9X)  C(*2  ,*2  .y  2 -y  2  *  ei)  de\d  didx\  dy'i  d*2dy' 2  j 


(4.  30) 

However,  the  plate  velocity  impulse  response  is  not  a  random  quantity;  hence 
the  ensemble  average  applies  only  to  the  turbulent  pressure  field.  Thus,  from 
(4.30)  and  (4.  29),  the  plate  velocity  cross  correlation  may  be  written  as 


<W*i-*2-yry2',r,2)-  111  II  I  Qpp  ,,2,y i ,y2 ,ci ~®i 2 ~ ^2^ 

foJoJojojo  JO 


.*'1  -y  1  .y  1 . 9 ,)  C(I  J  ,«  J  ,y  2  ,y  ' ,  ^2)  d  0,  d  dx'j  dy'  d* 2  dy \ 


(4.31) 
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From  Section  4. 1,  the  turbulent  boundary  layer  pressure  is  assumed  to 
be  a  homogeneous,  stationary  random  process.  Therefore,  the  turbulent 
pressure  cross  correlation  is  a  function  of  the  difference  between  the  spacial 
and  temporal  coordinates  rather  than  the  coordinates  themselves.  Hence, 

(4. 31)  may  be  written: 

mjojojojo  Jo 

<(* 2  •* 2  *y  j  <y j  •  e7)  d  6 1 d  e2  d*  i  dy \  *y \  • 

(4.32) 


<W,i-12-y,.y  rT) 


4.  2.  3  Plate  Velocity  Cross  Spectral  Density 

The  cross  spectral  density  is  defined  as  the  Fourier  Transform  of  the 
cross-correlation  function.  Hence,  the  plate  velocity  cross  spectral  density 
is  defined  as 


\[2n  J** 


Q^0(*,.*2.y,.y2.  T>«  dr. 


(4.  33) 


Multiplying  and  dividing  equation  (4.32)  by  and  substituting  into 

(4. 33)  gives  the  following  expression  for  the  plate  velocity  cross  spectral 
density : 


i 


.-*2-yi'y2^)-£££joYj[,{^j"Qpp(f>''T+0i-^)e‘‘“(T^,‘52,d(T  +  ^^)} 


<(*,.*[. y^y',.*,)*1^'  C(*2,Xj.yJ,y'2.(?2)e-1“92dfl1dflJdt'1dy/1d.jdy; 


(4.34) 

However,  the  term  in  brackets  is  easily  recognized  as  spp(«!f  .  7',  w);  there¬ 
fore,  (4.34)  may  be  rewritten: 
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W*  r*  i ’*1  «((  spp<£'. >!'.<»>  ("  C(* y ,, y',,  0,)  e““fl  1  d«, 

JoJoJoJo  Jo 

j"  <(*2,*'2.y2.y'2,e2)'luSi  dd:  d«',dy;dx'2dy'  . 

(1.35) 

Since  £(«,*',  y.y',  w)  is  zero  for  6  <  o  ,  the  semi-infinite  limits  in  i9(  and  02 
may  be  replaced  by  infinite  limits.  It  is  easily  shown  that 

(  i(*2,M2,y2y2,02)t‘iu6i  402  -  i£uH(«2,«'2,y2,y'Jt<j) 

J-oo 

and  that 

|  A6X  -  -iuH(s|ts'  ,yx,y\,-u)  . 

Hence, 

s**(li*,2’*|.yi' w)  "  m  wJS  pp  ( £  H(t ,  ,y ,  ,y',  ,-w) 

H(*2,«2,y2,y'  ,«)  d*' dy'd*2oy2 


(4.  36) 


(4.37) 


(4.  38) 

where  H  (»,*',  y ,  y <u)  is  as  defined  in  equation  (4. 16). 

The  expression  for  Spp(£',  r)  \  to)  was  presented  in  equation  (4.4).  Com¬ 
bination  of  (4.4),  (4.38),  and  (4. 16)  results  in 


no  no 


•<W1rI:->V>V“)  -  EL 


n(*1.y1)“q.(*2-y2) 
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Wv^-yryj^-^M-rH  £  £ 


2)  \  V'  ni  n'  I  7I  <j  ■  2  ’  2 


VJ  ::!  ::! 


h/*«rb/» 


>jujnjn 


f.o.  1 1  v-^i'cxlij  -  i|:>c-i(-  'c**2- '  cXIvj-vil* 


"mn^ry'i^q.^i-y^  d*')dy/1d*2dy'2^ 


>  1.2  V)  —  . 


(4.  39) 


where 


A  =0.-5  x  10*,o2f>(JU’S*. 


(4.40) 


Equation  (4.39)  may  be  rewritten  by  recalling  the  definition  of  amil(*.y)  from 
equation  (4.  8)  and  defining 

1.  =  »eK.;-.|'*sin^!l',in!lldlI'dl'  (4.41) 

",q  JoJo  “  “ 


and 


.flV-  l'c»'v2-y.1'  si 


nrryj  srry2 

sin - sin - dy'dy’  : 

b  b  1  : 


(4.42) 


nrry,  q»i;  sny. 


, . .  j  2  *  sin  - sin  — r—  sin - sin  ■  . — 
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It  remains  to  evaluate  I,  and  in  order  to  complete  the  solution 

mq  n » 

for  the  plate  velocity  case  spectral  density. 

From  (4.41),  I,  is  defined  as 

mq 


mq 


j^J*  e»0. 1 1  *<4Uch|  *2*  ‘  |l>e-l(w/lIcK 1 2  ■  » |  » 


'l  *2  .  , 

—  sin  — —  dx,  dx  ,  • 
a 


■l' 


(4.41) 


Note  that  the  Integrand  depends  on  the  absolute  value  of  «j-x'  ,  and  thus  the 
Integration  must  be  performed  over  limits  as  shown  in  Figure  6.  The  appro¬ 
priate  values  of  the  term  containing  the  absolute  value  for  each  area  of  inte¬ 
gration  are  also  presented  in  this  figure.  Thus,  by  utilizing  Figure  6,  equa¬ 
tion  (4.41)  becomes 


I,  *  dx j ^  dx'g 

+  f  e^/Ue><0  m  • 

(4.44) 

The  above  integral  may  easily  be  evaluated  from  standard  integral  tables,  and 
after  extensive,  but  routine,  simplification,  (4. 44)  becomes 

1.  -  — ^ —  is  1.0066  —  R  co»  (v  -  0.463  ») 

Ue  ■ 

r.  *  ((-il-t-l)"  -  ll  r„  iv„  „  .|v.i.,m4ff2  ,  * 

+  l 1  -  «,!  — i - T[R-e  q-R,e  -N2  —  co.(v  + 

•  hi-m 

_  mq,f  e-0.l»(»/Uc»[(_1)«  el(w./Uc  ♦  *,♦*'„>  +  (_1)qe-l«-  'Uc  ♦  v,+v„)]|  , 

(4.45) 

where  #- - - - j-j - - 

"■  -\0)  -  09B7(d7)] 4  °mn(^rj  •  <4-46» 
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(4.47) 

(4.48) 

o»(jf)' 

%  “  [(¥)’— (fjd' 

(4.49) 

5  is  the  Kronecker  delta. 

mq 

From  (4.42),  I,  is  defined: 

n  s 

n  b  b  d,'dy' 

(4.42) 

Again,  because  of  the  appearance  of  the  absolute  value  of  the  difference  of  the 
variables  of  integration  in  the  integrand,  the  integration  must  be  performed 
over  limits  as  shown  in  Figure  7.  By  applying  arguments  similar  to  those 
used  in  calculating  I ,  one  can  show  that 


4- 


—  [(-Dn(-i)‘ 


llli  -5  i+J!I!L[2  -|(-Dn 

n.  b  2 


(-l),|e°'7<“b/u«3[  • 


(4.50) 


Here  again,  5ns  is  the  Kronecker  delta.  It  shouldbe  noted  that  the  Kronecker 
delta  was  used  in  equations  (4.  45)  and  (4.  50)  to  conserve  space  rather  than 
writing  separate  expressions  for  m  =  q  and  m  ^  q  in  (4.45)  and,  similarly,  for 
n  =  s  and  n  4  s  in  (4.  50).  Hence,  in  order  that  no  confusion  arise  in  evaluating 
I,  and  l2  ,  the  effect  of  the  Kronecker  delta  should  be  considered  first. 

m  q  ns 
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Combining  (4.43),  (4.45),  and  (4.50),  one  finds 
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By  defining 
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equation  (4.51)  may  be  simplified  further.  The  resulting  expression  for  the 
plate  velocity  cross  spectral  density  is  given  by 
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(4.  60) 


4.  2. 4  Plate  Velocity  Spectral  Density 

The  power  spectrum  of  the  plate  velocity,  ,  y ,  to)  ,  may  be  defined 

in  terms  of  the  plate  velocity  cross  spectral  density  as  follows: 

<lV«i.y,.‘‘,>  -  (* | .* , ,y ,  ,y ,  ,a>)  .  (4*  61) 


Hence  from  (4.  60) 
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The  plate  velocity  power  spectral  density  should  be  a  real,  even  function. 
Hence,  upon  summation  equation  (4.  62)  must  be  real  in  order  that  it  be  a  valid 
solution  for  the  plate  velocity  spectral  density.  Substituting  for  V'm  from 

(4.59)  and  rearranging,  (4.62)  for  <a  <  1.256  ^  may  be  rewritten  as 

5 
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(4.63) 

Since  Gns  is  a  symmetrical  matrix,  that  is, 

G  -G  ,  (4.64) 

ns  sn 


the  first  summation  group  may  be  rewritten  as 
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C  (4. 65) 

By  applying  similar  arguments  to  the  last  three  summation  groups,  one  can 
show  that  (4.  62)  becomes 
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(4.  66) 
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for  oj  <  1.256 


U„ 


By  defining 


W  m  1.0066  —  5  R  cos  -  0.463")  cos  (A  -A  ) 

m  n  q  •  y  m  i]  in  in  inn  q  • 

U  c 


,,„>co,U  -A  >  K-l)-(-l)’  -  ll 

mnqi  m  q  2 

‘  —V  W 


2mq  n 

—  -  cos  . 

2  o  q 


[R  cos  (1/  +  A  -  A  )  -  R  cos  (i>  +  A  -  A  )) 

n>  q  m  n  q ■  q  m  q •  mn 


m<l'7  .0.  lIWi  t'c'l/  ,  \m  A*  .  i  \ 

- e  c  (-1)  cos  I —  +  p  +v  + A  -A  I 

y  1 IJ  o  ■»  mn  1*1 

(-l)q  ros/—  +  v  +  i  ♦  A  -  A  \"1 
1  [J  q  m  q*  mnl  I 


(4.  67) 


and  noting  that  arguments  applied  to  equation  (4.  65)  for  m  1.256  are 

l > 

U„ 

also  valid  for  the  frequency  range  t*.  •  1.256  —  ,  one  may  write  the  plate 

8 

velocity  power  spectral  density  as 
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-I.  :*  Development  of  the  Cavity  Acoustic  Pressure  Cross  Spectra)  Density 
4.  9.  1  Cavity  Acoustic  Pressure  due  to  an  Arbitrary  Plate  Velocity  Distribution 
The  equations  governing  acoustic  phenomena  are  the  momentum  equation, 


4- 


vp<  -  o, 


(4.  69) 


the  continuity  equation, 


dp 

- '  +  p  V  .  U  m  0  , 

,h  P*0 


(4.70) 


and  the  equation  of  state, 

Pt«c 2  Pt.  (4.71 

It  should  be  noted  that  the  subscript  "zero"  in  the  above  equations  refer 
to  time  average  quantities,  whereas  non-subscripted  quantites  are  the  instan¬ 
taneous  or  fluctuating  components. 

Reference  to  Figure  1,  for  an  arbitrary  plate  velocity  $(i,y  ,t) ,  indi¬ 
cates  that  equations  (4.  69),  (4.  70),  and  (4.  71)  must  be  solved  subject  to  the 
following  boundary'  conditions: 


11^(0,  y,z,  t)  *0  (•) 

uj(a,  y,  z,  t)  -0  (b) 

u  («,  0,  z,  0  =0  (c) 

y  (4.72) 

uy(i,  b,  z,  t)  "0  (d) 

uj(z,y,0,t)  -0  (e) 

u^z.y.-d,  l)  -  <£(*,y,t)  (f)  ■ 


Equation  (4.  69)  was  derived  on  the  basis  of  an  inviscid  fluid.  Hence,  the 
acoustic  field  may  be  assumed  irrotational,  and  the  acoustic  phase  velocity 
may  be  defined  in  terms  of  a  velocity  potential  U'<)  as  follows 


u(l,  y,  Z,  t)  -  V  0(i,  y,  z,  t)  . 


(4.79) 
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Equation  (4.73)  only  specifies  o  to  within  an  arbitrary  function  of  time. 
Hence,  in  order  that  <!>  be  uniquely  specified,  it  is  necessary  to  further 
define 

d4<x,y,  z.  0 

Pt(x,y,z,t)  - - - .  (4.74) 


Substitution  of  (4.  73)  and  (4.  74)  into  (4.  69)  identically  satisfies  (4.  69).  It  re¬ 
mains  only  to  satisfy  (4.  70)  and  (4.  71).  Substituting  (4.  71)  into  (4.  70),  one 
finds  that 


1  *P. 

ci  dt 


•  U 


0. 


(4.75) 


Substitution  of  (4.  74)  and  (4.  73)  into  (4.  75)  results  in  the  scalar  wave  equa¬ 
tion  in  0:  that  is, 


V2  0  - 


1 

c2  dt1 


-0. 


(4.  76) 


Separation  of  variables  techniques  yields  a  solution  to  (4.  76): 


0(i,y,  z,  0  -  [C(  sin  k>»  +  Cj  cos  k^]  [C}  sin  kyy  +  C4  cos  kyyl 

(C,  sin  k  z  +  C  cos  k  z]  C,  e,l,c,  , 

'  to  Z  7 


where 


k2  +k2  +  k2  -k2. 

i  y  i 


(4.  77) 


(4.  78) 


Application  of  the  boundary  conditions  (4.  72a)  through  (4.  72e)  to 
(4.77),  by  means  of  (4.73),  yields 


equation 


(4.  79) 
(4.  80) 
(4.81) 


where  (  and  t  are  integers.  Further,  by  defining 


.'ifi 


«<)  »  kc 


(4.  82) 


and  by  utilizing  (4.  79),  (4.80),  and  (4.81),  it  is  seen  that 

0(«,y,  t,  t)  -  V  C2  cosilic  cos— ^Cficos(L  z)C  e1  (4.89) 

j»o  i  •  I  b  if 

t-o 

Combining  all  the  constants,  one  may  assume  that  the  solution  for  tM* ,  y .  z ,  i) 


takes  the  form 
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i*0  *  b  JU,  '«  \2rr 

l*o 


Substitution  of  (4.83a)  into  (4.76)  readily  shows  that  (4.83a)  satisfies  the 


wave  equation,  and  by  applying  the  above  arguments,  boundary  conditions 


(4.  72a)  through  (4.  72e)  are  satisfied.  It  therefore  remains  only  to  satisfy  the 


final  boundary  condition  (4.72f):  that  is, 

0(l,  y,  t)  **  cos  cos — -  f  Y..(ai)|-k  sin  Ik  (-d)l(ell*)t  -M-  .  (4.84) 
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Let 


X  ((j)=k  (oj)  Y  ,(<.))  sin(k  d) 
il  *  j*  11  'it 

Equation  (4.  84)  may  be  rewritten  as 


(4.  85) 


,,  .  v"  i nx  uy  f  v  i  \ 

0(l,y,t)-  >  cos -  cos -  I  X  (oj)  e  - . 
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(4.  86) 


Equation  (4.  86)  is  essentially  a  two-dimensional  Fourier  cosine  series  in  *  and 


y  ;  hence  X(<,>)  may  be  evaluated  by  determination  of  the  Fourier  coefficients 


of  (4.  86).  By  utilization  of  the  orthogonality  principle,  it  may  be  shown  that 
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_J _ 
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(4.  87) 
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Transformation  of  (4.  87)  yields 

X  M  * - - -  C  C  y, «)  cot  - — cos— —e'iU1'  dx  dy— (4.88) 
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Thus,  from  (4.  85) 
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and  using  this  result  in  (4.  83a)  yields 
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(4. 90) 


The  cavity  acoustic  pressure  for  an  arbitrary,  deterministic  plate 


velocity  is,  from  (4.74)  and  (4.90), 
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(4.91) 


4.3.  2  Cavity  Acoustic  Pressure  Cross  Correlation 


The  cavity  acoustic  pressure  cross  correlation  is  defined  as 


?i.(ll'l2')rl-)r2-,l-l2>‘i  «2)  "  ED,a,*pyrzl-,l)P.(*2*>'2'*2-,2)]  . 


(4.92) 


From  (4.  SI)  and  (4.92),  by  utilizing  the  lav  vledge  that  the  only  random  quan¬ 


tity  in  the  expression  for  pa(*  , y .  z ,  t)  is  the  plate  velocity,  one  may  write  the 


acoustic  pressure  cross  correlation  as 
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(4.  93) 


From  (4.28)  and  (4.32),  however, 

F.  [«5(*  i  .y  i  .«i>  ^<*,.y,.«2>]  =  Q*,  <*,.*2.y,.y2.T) 

Therefore,  (4. 93)  becomes 
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(4.  95) 


It  is  now  of  interest  to  examine  the  following  portion  of  (4.  95): 


By  recalling  that  t  c ,  -  c  ]  .  one  may  rewrite  (4.  96) : 


(4.  96) 


(4.97) 
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However,  from  (4.33) 
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Therefore,  (4.  97)  becomes 
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Finally,  since 

|  e‘,,‘  -  2f7S(oj  ♦  H) 

equation  (4.98)  becomes 

(2,r)’ 1 5W*r*2>>,i*y2*fl>  «<<*>  4  n>. 

Hence,  (4.  93)  becomes 
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(4. 100) 
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(4. 101) 


Integration  over  <u  .  and  use  of  the  identity 


t  ,  ~  t  I 


lead  to  the  conclusion  that  equation  (4. 101)  may  be  written: 
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(4.  102) 


It  should  be  noted  that,  from  (4.  78),  (4.  80),  (4.  81),  and  (4.82), 


k, 


to) 

! 


(4. 103) 


Thus, 


k,  (<■)  k,  (4.104) 


and  the  cavity  acoustic  pressure  cross  correlation  may  be  written: 
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4.3.3  Cavity  Acoustic  Pressure  Cross  Spectral  Density 

The  cavity  acoustic  pressure  cross  spectral  density  is  defined  as  the 
Fourier  transform  of  the  cross  correlation:  that  is, 
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Thus,  from  (4. 105) 


Hr,2  '*  ' 


jrrij  tnyl  irry? 


•  O  'I  '1  2  7  2 

saJ(xi  -y^**  i  j  '•  >  — 7-7  L  2^  cos  — co* -r-co*  —  to*-7— 


rr a'b‘  i  =  <>  f  =  <> 


(££ 


!  =  n  , :  o 

13  ‘  COS  If,  7.  cos  k,  7 
_  It  1  ft  1 


d  *  ^0  )(1  •  <^nf)(l  •  >>nf)(l  *  ^0,)k,|((n)k/(|(i7)  sin  k,  d  s,n  k, f  d 


[ffff 


j  it  * ,  Iffy,  rffx,  iffy : 

Srf,  Is  |  ,x  ,  ,y  |  ,y  ,  .11)  cos - cos  - cos  - -  cos - 

’  *  i  i  i  .  a  «  a  L. 


dx  jdyjdx^iy  ,1 


«-,li7dil  c“  7  dri. 
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If  the  expression  for  the  plate  velocity  cross  spectral  density  from  (4.  60)  is 
now  substituted  into  (4. 107)  and  subscripts  are  reassigned  such  that  j,  k,  m  , 

and  n  apply  to  the  acoustics  problem  and  q,  r,  s,  and  i  apply  to  the  plate, 
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From  the  above  definitions,  (4. 108)  may  be  rewritten  as 
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The  form  of  (4. 114)  indicates  that  the  integration  over  il  and  the  space  vari¬ 
ables  x2,  y, ,  and  y2  may  be  performed  separately.  This  is  easily 
seen  by  rearranging  (4. 114)  as  follows: 
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Since  the  inverse  Fourier  transform  is  detined  as 
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the  integration  over  frequency  of  the  bracketed  terms  in  (4. 115)  results  in 
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When  (4.110)  li  iniiorted  Into  (4.  US)  and  the  integration  over  7  is  performed, 
(4. 115)  becomes 
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It  remains  to  perform  the  integration  over  the  spaclal  coordinates.  Using 
(4. 109)  and  (4. 112),  by  means  of  standard  Integration  techniques,  one  findB 
that 
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From  (1. 117)  and  (4. 119),  the  cavity  acoustic  pressure  cross  spectral  density 
is 
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(4. 120) 

where  Gjknin<«l,.J,yI,y2>,  *,«„(*,. «V*).(,,,>>.XU)  and  Kjkmnqfi,  are  defined 
in  equations  (4. 112),  (4.113),  (4.110),  (4.111),  and  (4. 119),  respectively. 
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4.3.4  Cavity  Acoustic  Pressure  Spectral  Density 

The  cavity  acoustic  pressure  spectral  density  is  defined  as 
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Hence,  from  (4. 120), 
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Hence,  using  (4. 110),  (4.  Ill),  (4. 123),  and  (4. 124)  in  (4. 122),  and  dropping 
the  subscripts  on  the  spacial  coordinates,  one  may  express  the  acoustic  spec¬ 


tral  density  as 
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(4. 126) 

As  In  the  case  of  the  plate  velocity  spectral  density,  the  cavity  acoustic 
pressure  spectral  density  must  be  a  real,  even  function  of  frequency.  Thus, 
by  substituting  v  ft|  from  (4. 59)  and  by  investigation  of  the  summation  in  the 
frequency  range  u  $  1.256  ^  .  arguments  similar  to  those  applied  in  Section 
4. 2. 4  along  with  equations  (4. 123),  (4. 124),  and  (4. 125)  show  the  first  sum¬ 
mation  group  to  be 
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Hence,  tho  first  summation  group  is 
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Similar  arguments  applied  to  the  remaining  summation  groups  lead  to 
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4.4  Non-DimenslonpJ  Spectra 

The  results  of  Sections  4.  2.4  and  4.3.4  are  more  useful  from  a  practical 
viewpoint  if  they  are  expressed  in  a  non-dimensional  form.  The  following  sec¬ 
tions  are  devoted  to  expressing  the  plate  velocity  and  cavity  acoustic  pressure 
spectral  densities  in  non-dimensional  forms. 

4.4. 1  Dimensionless  Plate  Velocity  Spectral  Density 

The  plate  velocity  spectral  density  (equation  (4. 68))  is  a  function  of  the 
following  parameters  and  dimensions: 
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Note  that  only  three  dimensions  (force,  length,  and  time)  are  required  to  spe¬ 
cify  the  eleven  Input  parameters.  It  is  therefore  necessary  to  Belect  a  length, 
force,  and  time  scale  characteristic  of  the  system.  Although  there  are  many 
combinations  of  the  eleven  input  parameters  which  result  in  dimensions  of 
force,  length,  or  time,  the  following  were  selected  as  characteristic  dimen¬ 
sions  because  it  was  felt  that  they  had  the  strongest  effect  on  the  plate  velocity 
spectral  density: 


a  -  streamwise  plate  dimension  (L) 

—  -  time  for  an  eddy  to  traverse  the  plate 

uc  in  the  streamwise  direction  (T) 
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fictitious  force  based  on  the  mass  of  the 
plate  (F) 
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By  utilizing  those  characteristic  dimensions,  It  is  possible  to  specify 
dimensionless  input  parameters  as  follows: 
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Further,  the  spacial  and  frequency  variables  may  by  similarly  non- 
dimensionalized  as  follows: 


(4. 130) 


Based  on  the  definition  of  (4. 129)  and  (4. 130),  the  quantities  defined  in  equa¬ 
tions  (4.46)  through  (4.49),  (4.  52)  through  (4.  58),  and  equation  (4.67)  may  be 
rewritten  in  dimensionless  terms  as  follows: 
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Further,  by  defining 
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and  utilizing  the  above  definitions  and  equation  (4.  68),  the  non-dimensional 
form  of  the  plate  velocity  spectral  density  can  easily  be  shown  to  be 
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4. 4. 2  Dimensionless  Cavity  Acoustic  Pressure  Spectral  Density 

The  cavity  acoustic  pressure  spectral  density  is  a  function  of  all  of  the 
parameters  affecting  the  plate  velocity  spectral  density  (see  Section  4. 4. 1) 
plus  the  following  additions: 
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In  addition,  there  is  one  more  variable  in  the  case  of  the  acoustic  spectral 
density,  namely,  the  length  dimension  (2)  perpendicular  to  the  lane  cl  !Vo 
plate. 

Since  the  plate  vibration  provides  the  excitation  for  the  acoustic  pressure 
in  the  cavity,  it  would  seem  reasonable  to  select  the  same  characteristic 
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dimensions  for  non-dimonslonallzatlon  of  the  cavity  acoustic  pressure  spectral 
density  as  were  used  for  the  plate  velocity  spectral  density.  The  dimension¬ 
less  input  parameters  for  the  cavity  acoustic  pressure  spectral  density  are, 
therefore,  those  presented  in  equation  (4. 129)  plus  the  following  additions: 
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Likewise,  the  variables  are  those  presented  in  equation  (4. 130)  plus 

(4.141) 


Since  the  non-dlmensionalization  has  been  performed  with  the  same  character¬ 
istic  dimensions  as  in  the  case  of  the  plate  velocity  spectral  density,  the 
dimensionless  quantities  defined  in  equations  (4. 131)  to  (4. 137)  are  also  appli¬ 
cable  to  the  cavity  acoustic  pressure  spectral  density.  Noting  that  Kjkmnqtsi 
is  already  dimensionless,  one  needs  oiJy  to  define 


-|‘- 


<V)k!  («♦) 


sin  k*  d  sin  k +  d+  , 
*ik  zmn 


(4. 143) 


and 


=  - - — - - - 

Qyuc5 


(4. 144) 

in  order  to  specify  the  dimensionless  form  of  the  cavity  acoustic  pressure 


spectral  density.  By  use  of  equations  (4. 128),  (4. 131)  through  (4. 137),  and 
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(4. 142)  through  (4. 144),  the  dimensionless  form  of  the  cavity  acoustic  spec 
tral  density  can  easily  be  bhown  from  (4. 128)  to  be 


TJr  (w+)Pr  P«+  <«*>  R  *(w+)R;  U>*) 


u*8'  <  1.932  , 


a>*8*  >  1.932  . 


(4. 145) 


V 


RESULTS 

The  following  section  presents  plate  velocity  spectra  and  cavity  acoustic 
spectra  calculated  from  equations  (4. 139)  and  (4. 145)  by  means  of  a  digital 
computer  for  selected  data  cases.  The  purpose  of  these  calculations  was 
threefold.  First,  the  calculated  plate  velocity  spectral  density  may  be  com¬ 
pared  with  existing  measurements  of  plate  response  to  turbulent  excitation, 
thereby  providing  a  check  on  the  validity  of  the  analytic  solution  for  the  plate 
velocity  statistics.  Secondly,  the  effects  of  all  major  parameters  are  not 
readily  discernible  from  equations  <4. 139)  and  (4. 145),  and  it  was  felt  that 
calculated  spectra  for  selected  cases  would  better  define  these  effects.  Fi¬ 
nally,  it  was  necessary  to  check  the  range  of  validity  of  the  assumption  that 
the  cavity  acoustic  pressure  had  negligible  effect  in  the  excitation  of  the 
plate. 

The  spectral  densities  were  selected  for  computation  over  the  cross 
spectral  densities  because  they  were  most  easily  computed  and  gave  the  de¬ 
sired  results:  namely,  the  effect  of  major  parameters  on  the  cavity  acoustic 
pressure  statistics,  and  the  range  of  validity  of  the  above  stated  assumption. 

The  computations  were  performed  partially  on  the  IBM  704  digital  com¬ 
puter  at  the  U.  S.  Navy  Underwater  Sound  Laboratory,  and  the  remainder  on 
the  IBM  7090  digital  computer  at  the  David  Taylor  Model  Basin. 
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5. 1  Results  of  Plate  Velocity  Spectral  Density  Computations 

The  plate  velocity  spectral  density  was  computed  from  the  dimensionless 
equation  (4. 139).  Computations  were  performed  on  the  IBM  704  digital  com¬ 
puter  by  an  iterative  process.  In  all  cases,  the  summation  terms  in  equation 
(4. 139)  appeared  to  converge  monotonically,  and  the  solution  was  accepted  if 
two  successive  iterations  agreed  within  three  percent.  ThiB  tolerance  was 
selected  as  a  compromise  between  desired  accuracy  and  reasonable  computa¬ 
tion  time. 

The  dimensionless  plate  velocity  spectrum  was  generated  by  computing 
the  spectral  density  at  a  dimensionless  frequency  one  cycle  above  each  dimen¬ 
sionless  plate  natural  frequency,  and  at  three  points,  equally  spaced  in  fre¬ 
quency,  between  each  pair  of  adjacent  natural  frequencies.  Peculiarities  of 
the  particular  computer  program  used  for  the  generation  of  the  frequencies  at 
which  the  spectral  density  was  to  be  computed  were  responsible  for  the  spec¬ 
tral  density  being  computed  at  a  dimensionless  frequency  one  cycle  higher  than 
the  plate  natural  frequency.  The  resulting  error  in  the  computed  plate  veloc¬ 
ity  spectrum  at  plate  resonances  should  be  very  small,  except  at  very  low 
frequencies  or  for  very  small  values  of  damping. 

5. 1. 1  Comparison  with  Existing  Experimental  Data 

In  order  to  check  the  validity  of  the  analytically  derived  plate  spectral 
density,  it  was  desirable  to  compare  the  computed  plate  velocity  spectral  den¬ 
sity  from  equation  <4. 139)  to  available  experimental  data.  Unfortunately,  all 
existing  experimental  information  deals  with  fixed  edge  plates  whereas  equation 
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(4. 139)  pertains  only  to  a  simply  supported  plate.  However,  the  spectral 
levels  for  a  simply  supported  plate  should  be  of  the  some  order  of  magnitude 
as  those  for  a  plate  with  clamped  edges,  even  though  the  plate  natural  frequen¬ 
cies  will  be  different  in  the  two  cases.  It  was  therefore  decided  to  compare 
the  general  level  of  the  spectrum  predicted  by  equation  (4. 139)  with  available 
data  for  a  clamped  edge  plate. 

Some  of  the  tro^t  pypriWmont'*!  dnt?  on  the  flmv  excita¬ 

tion  of  flat  plates  are  those  of  Bull  rial.  [15],  Bull  measured  the  displacement 
spectral  density  of  3. 5  *  3. 5  inch  steel  plates  of  varying  thicknesses  to  turbu¬ 
lent  boundary  layer  excitation  in  a  wind  tunnel.  Since  no  acoustic  cavity  was 
placed  behind  the  plate,  the  effects  of  acoustic  back  pressure  are  negligible. 
Figure  8  presents  Bull's  measurements  of  the  response  of  a  3. 5  x  3. 5  *  0. 010 
inch  plate  to  air  flowing  at  a  free  stream  velocity  of  539  ft/sec.  The  meas¬ 
ured  displacement  thickness  was  0. 172  inch.  The  measurement  was  taken  at 
a  position  0. 95  inch  upstream  from  the  trailing  edge  of  the  plate  and  centered 
in  the  lateral  direction.  Bull  measured  the  modal  damping  for  the  first  two 
plate  modes  and  found  them  to  be  3.  8  percent  critical  and  0. 5  percent  criti¬ 
cal,  respectively.  He  states,  however,  that  the  measurements  were  accurate 
only  to  +23  percent. 

Figure  9  presents  the  computer  solution  to  equation  (4. 139),  converted 
to  the  ratio  of  displacement  spectral  density  to  turbulent  pressure  spectral 
density,  for  input  data  corresponding  to  Bull's  experiment.  The  dimension¬ 


less  input  parameters  were 
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M  -  0.099 
11  *  -  1.54  (■■turned) 

5*  -  0.0491 
b*  -  1.0 
D*  -  0.001745 

r  ♦ 

—  «  0.005  (assumed)  . 
r  * 

an 

Since  g/ b  and  r*  are  functions  of  D*  and  b*  they  need  not  be  independ¬ 
ently  specified.  It  should  also  be  noted  that  u\  which  specifies  the  ratio  of 
the  free  stream  velocity  to  the  convection  velocity,  is  assumed  to  be  constant. 
This  assumption  is  not  true  in  fact,  since  the  value  of  U4  is  a  function  of 
gj5*/U0  as  shown  in  Figure  4.  However,  computations  are  greatly  simplified 
by  assuming  u*  to  be  constant.  Finally,  note  that  the  modal  damping  is 
assumed  to  be  constant  with  frequency  and  equal  to  the  value  measured  by  Bull 
for  the  2-1  mode.  In  general,  the  modal  damping  is  not  constant,  but  this 
assumption  eased  computational  difficulties.  In  the  computation  of  the  spec¬ 
trum  of  Figure  9,  the  spectral  densities  were  computed  at  radial  frequencies 
0. 005  cycle  (that  is,  equal  to  the  damping  -  critical  damping  ratio)  above  the 
plate  natural  frequencies,  and  at  three  points  equally  spaced  in  frequency 
between  each  adjacent  paifr  of  natural  frequencies.  This  procedure  resulted 
in  more  accurate  spectral  levels  for  the  low  values  of  damping  being  used. 

Figure  10  compares  the  theoretically  derived  spectrum  for  the  simply 
supported  edge  condition  (Figure  9)  with  Bull's  measurements  for  the  fixed 
edge  condition  (Figure  8).  Note  that,  as  expected,  the  natural  frequencies  are 
not  in  agreement  because  of  the  different  edge  conditions  imposed  on  the  two 
cases.  Hence,  the  spectra  for  the  two  cases  are  not  in  exact  agreement. 
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However,  the  envelopes  of  both  spectra  are  In  substantial  agreement,  especi¬ 
ally  at  the  higher  frequencies.  This  fact  lends  credence  to  this  analynis  and 
to  the  ability  of  equation  (4. 139)  to  predict  correctly  the  velocity  spectral 
density  of  simply  supported  plates. 

For  the  higher  plate  modes,  it  seems  reasonable  that  the  modal  stiffness 
and  mass  become  less  and  less  affected  by  the  plate  edge  conditions.  Hence, 
for  this  linear  system,  it  would  further  seem  reasonable  that  a  line  connecting 
the  major  spectral  peaks  in  the  clamped  edge  case  should  be  very  similar  in 
shape  and  amplitude  to  one  connecting  the  major  spectral  peaks  in  the  simply 
supported  case,  especially  at  the  higher  frequencies.  For  purposes  of  discus¬ 
sion  here,  such  a  line  will  be  defined  as  a  "peak  spectrum. "  Figure  11  pre¬ 
sents  the  "peak  spectra"  derived  from  Figures  8  and  9.  The  agreement  in 
shape  is  excellent  above  700  cps,  and  the  amplitudes1  agree  within  5  db.  This 
discrepancy  in  amplitude  may  easily  be  due  to  discrepancies  between  the 
assumed  value  of  the  damping  used  in  the  computed  spectrum  and  the  damping 
present  at  higher  modes  in  Bull's  experiment. 

From  the  above,  it  appears  that  the  plate  velocity  spectral  density  pre¬ 
dicted  by  equation  (4. 139)  is  valid  for  the  simply  supported  plate  excited  by 
turbulent  boundary  layer  pressure  fluctuations  for  cases  in  which  the  acoustic 
environment  has  negligible  effect.  It  may  further  be  concluded  that  the  "peak 
spectrum"  of  a  flow  excited  plate  is  essentially  independent  of  the  plate 
boundary  conditions  at  frequencies  above  the  first  few  plate  natural  modes. 
From  an  engineering  standpoint,  this  conclusion  allows  an  estimation  of  the 
maximum  flow  excited  plate  vibration  level  .  bove  a  certain  frequency  from 
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equation  (4. 139}  without  regard  to  boundary  conditions.  Determination  of  this 
lower  frequency  limit  will  be  treated  in  the  next  section. 

5. 1.  2  Effect  of  Major  Parameters  on  Plate  Velocity  Spectrum 

As  stated  in  Section  II,  the  objective  of  this  study  is  to  provide  the  sub¬ 
marine  sonar  systems  designer  with  information  concerning  the  acoustic 
environment  of  sonar  transducers  and  the  major  parameters  which  may  affect 
this  environment.  Therefore,  the  data  cases  used  to  show  the  effects  of  the 
major  parameters  on  the  plate  velocity  spectrum  were  selected  with  the  sub¬ 
marine  application  in  mind.  It  must  be  stressed  that  in  all  cases  the  effects 
of  any  plate-generated  acoustic  pressure  have  been  neglected. 

The  major  input  parameters  for  the  dimensionless  plate  spectral  density 
are  presented  in  equation  (4. 129).  For  the  purposes  of  this  discussion,  the 
variables  stated  in  equation  (4. 130)  may  also  be  considered  major  parameters. 
From  the  definitions  of  <«/  and  r*  ,  it  can  be  shown  that  these  parameters 
are  functions  only  of  D+  and  b*.  Further,  U*  is  assumed  to  be  equal  to  1.54 
throughout  this  analysis.  Hence,  the  only  parameters  which  need  be  specified 
in  order  to  obtain  a  solution  to  (4. 139)  are  m,  5+,b+,r+,D+,i+,y  +  ,  and  the  range 
of  interest  of  u*  . 

By  inspection  of  equation  (4. 139)  and  equations  (4. 131)  through  (4. 137),  it 
is  evident  that  effects  of  b\r*,D  +  ,x+,y*,  and  o»*  are  obscured  by  the  summa¬ 
tion  term  in  (4. 139).  However,  the  parameters  M  and  8*  occur  only  outside 
of  the  summation,  so  that  the  effect  of  these  parameters  on  the  dimensionless 
plate  velocity  spectral  density  are  immediately  obvious.  For  most  submarine 
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applications,  b*  will  vary  between  0.5  and  1.0,  and  it  was  felt  that  varia¬ 
tions  in  b*  over  this  range  would  change  the  plate  natural  frequencies  but 
would  not  significantly  alter  the  general  level  of  the  spectrum.  Hence,  the 
effect  of  plate  shape  (b*)  was  not  investigated  in  this  study.  This  omission 
leaves  only  the  damping  (r*) ,  the  plate  rigidity  (D*) ,  plate  point  of  interest 
(x  +  ,y+) ,  and  frequency  (&/)  as  parameters  to  be  investigated.  The  effects  of 
each  of  these  parameters  will  be  discussed  below. 

5. 1. 2. 1  Effect  of  Plate  Coordinates 

Dimensionless  input  parameters  which  fall  in  the  range  of  interest  for 

submarine  applications  are  as  follows: 

M  =  18.9 
D+  =  10.3 
S*  =  2.31  x  10' 2 
b+  =  0.6667 

r+  =  0.1  r*  (assumed), 
m  n 

The  position  of  maximum  plate  vibration  is  usually  of  greatest  interest  as  far 
as  plate  coordinates  are  concerned.  It  seems  reasonable  that  the  maximum 
level  would  occur  at  the  center  of  the  plate.  Hence,  the  effect  of  plate  coordi¬ 
nates  was  investigated  over  the  range  of  plate  coordinates  0.2  <x+  <0.5,  y  +  <0.5  . 
The  frequency  range  of  interest  was  selected  to  be  10  <  <  3000. 

Figures  12,  13,  and  14  present  dimensionless  plate  velocity  spectra  for 
the  above  stated  case  with  plate  coordinates  (*  +  ,y+)  of  (1/5,  1/3) ,  (1/3, 1/3) , 
and  (1/2, 1/3) ,  respectively.  These  three  spectra  show  significant  differ¬ 
ences  in  shape  below  a  dimensionless  frequency  of  1000  because  of  the  points 
oi  measurement  falling  on  different  plate  modal  nodes  or  antinodes.  A  table 


of  the  plate  dimensionless  natural  frequencies  below  1000  dimensionless 
cycles  is  presented  in  Table  1.  Note,  for  example,  that  **  =  1/2,  y*  =  1/3 


TABLE  1 

PLATE  DIMENSIONLESS  NATURAL 
FREQUENCIES  BELOW  u*  =  1000 


M 

D* 

Input  Parameters 

-  18.9  5*  =  0.0231 

=  10.3  b*  =  0.6667 

r*  =  0.1  r* 

ran 

(o 

mn 

m 

i 

2 

3 

i 

102.9 

316.8 

673.1 

2 

198.0 

411.8 

768.1 

3 

356.3 

570.2 

926.5 

4 

5 

578.1 

863  1 

791.9 

places  the  point  of  interest  at  the  plate  center;  hence  one  would  not  expect  to 
observe  u{2,  u,  <i>22,  etc.  However,  this  position  is  antinode  for 
uj ,,  ojjj.  wj,,  etc.  Figure  14  clearly  shows  resonances  at  the  frequen¬ 
cies  (from  Table  1)  associated  with  the  1-1,  1-3,  3-1,  and  5-1  modes, 
whereas  the  1-2,  2-1,  4-1,  and  4-2  modes  are  not  discernible.  Similar 
effects  may  be  noted  in  Figures  12  and  13.  This  fact  further  supports  the 
validity  of  equation  (4. 139). 

Consider  now  the  "peak  spectra"  derived  from  Figures  12,  13,  and  14 
above  the  first  plate  natural  frequency.  A  comparison  of  these  "peak  spectra" 
is  presented  in  Figure  15.  Inspection  of  Figure  15  shows  negligible  effect  of 
plate  coordinates  on  the  dimensionless  plate  velocity  "peak  spectra  density" 
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above  the  first  plate  mode.  It  Is  concluded,  therefore,  that  plate  coordinates 
in  the  range  0.2  <  <  0.3,y*  «=  0.333  have  no  effect  on  the  dimensionless  plate 

velocity  "peak  spectral  density. "  This  conclusion,  In  addition  to  the  observa¬ 
tions  made  above  concerning  Figures  12,  13,  and  14,  supports  the  further 
conclusion  that,  above  the  first  natural  frequency,  the  effect  of  changing  the 
point  of  observation  on  the  plate  (plate  coordinates)  is  a  change  in  shape  of  the 
dimensionless  plate  velocity  spectrum  without  changing  the  general  level  of 
the  spectrum. 

5. 1.  2.  2  Effect  of  Plate  Damping 

By  utilizing  the  same  input  parameters  listed  in  Section  5. 1. 2. 1,  the 
dimensionless  plate  velocity  spectrum  was  computed  for  two  additional  values 
of  damping  at  the  midpoint  of  the  plate  (**  =  1/2 ,  y*  =  1/3)  over  the  frequency 
range  10  <  <*/  <  3000  .  These  additional  values  of  damping  were  0.01  t*  and 

"  =  -  cmn 

0. 05  r*  The  resulting  dimensionless  plate  velocity  spectra  for  these  cases 
are  presented  in  Figures  16  and  17,  respectively.  Thus,  Figures  14,  16,  and 
17  show  the  effect  of  damping  alone  on  the  dimensionless  plate  velocity  spec¬ 
tral  density.  Comparison  of  these  figures  shows  that  an  increase  in  damping 
decreases  the  spectral  density  in  the  region  of  the  natural  frequencies  present 
in  the  spectra,  as  would  be  expected.  The  amount  of  this  reduction  in  spectral 
level  is  most  easily  investigated  by  inspection  of  the  dimensionless  "peak 
spectra"  constructed  from  Figures  14,  16,  and  17.  A  comparison  of  the 
"peak  spectra"  for  the  three  damping  cases  is  presented  in  Figure  18.  It 
should  be  noted  here  that  the  first  plate  dimensionless  natural  frequency 
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1  932 

falls  above  to*  - - ,  and  that  the  "peak  spectra"  presented  in  Figures  15  and 

5 

18  only  apply  above  ojJ, .  Hence,  any  conclusion  drawn  from  the  dimensionless 
"peak  spectra"  apply  only  to  dimensionless  frequencies  above  1.932/5*. 

Figure  18  shows  the  effect  of  damping  on  the  dimensionless  plate  velocity 
spectral  density  at  the  plate  natural  frequencies.  Note  that  the  slope  on  the 
"peak  spectrum"  is  also  slightly  changed  as  the  damping  is  varied.  This 
effect  may  be  a  result  of  the  computational  procedure  of  computing  the  dimen¬ 
sionless  spectral  density  at  a  dimensionless  frequency  one  cycle  higher  than 
the  plate  natural  dimensionless  frequency.  Since  the  slope  change,  as  noted 
above,  was  slight,  it  was  not  felt  that  this  effect  warranted  a  complete  investi¬ 
gation.  The  effect  of  damping  on  the  slope  of  the  "peak  spectrum"  will  be 
discussed  as  it  appears  in  Figure  18.  However,  it  should  be  kept  in  mind  that 
this  effect  may  have  no  physical  basij. 

Figure  19  presents  the  effect  of  damping  on  the  "peak  spectrum"  over  the 
range  of  damping  o.oi  < -I—  <0,1  for  the  dimensionless  frequency  of  200, 

r  * 

cm  n 

which  is  equivalent  to  W*S*  =  4.62 .  This  effect  may  be  expressed  as 


(5.  i) 

mn  / 

It  is  obvious  from  (5. 1)  that  the  "peak  spectrum, "  and  hence,  the  dimension¬ 
less  plate  velocity  spectral  density  at  the  plate  natural  frequencies,  varies 
inversely  as  the  second  power  of  plate  damping  at  the  dimensionless  frequency 
of  4.62/6* 

Figure  20  presents  the  variation  in  slope  of  the  dimensionless  "peak 
spectrum*'  with  damping  over  the  same  range  of  damping.  This  effect  may 
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be  described  by 


From  (5.  2),  it  is  seen  thpi,  for  dimensionless  frequencies  above  the  first 
plate  dimensionless  natural  frequency  and  above  1.932/5* ,  the  slope  of  the 
dimensionless  "peak  velocity  spectrum"  increases  with  the  one-half  power  of 
damping. 

From  the  above,  it  may  be  concluded  that  the  plate  damping  has  a  pro¬ 
nounced  effect  on  the  dimensionless  plate  spectral  density  at  dimensionless 
frequencies  corresponding  to  the  plate  natural  frequencies  present  in  the  spec¬ 
trum.  The  major  effect  of  damping  is  to  decrease  the  dimensionless  spectral 
density  at  these  frequencies  according  to  approximately  the  inverse  second 
power  of  the  damping.  This  conclusion  was  also  reached  by  Bull  et_al.  [l5]  in 
their  investigation  on  the  flow  excited  vibrations  of  plates.  This  lends  further 
support  to  the  validity  of  this  analysis. 

5. 1.  2. 3  Effect  of  Plate  Rigidity 

By  use  of  the  values  of  M ,  5* ,  b* ,  and  r*  prescribed  in  Section  5. 1. 2. 1, 
the  dimensionless  plate  velocity  spectral  density  was  computed  at  the  plate 
center  (x*  =  i/2,  y*  =  1/3)  over  the  dimensionless  frequency  range  of  11  to 
3000  for  the  additional  values  of  D*  of  5, 15  and  15.45.  The  results  of  these 
computations  are  presented  in  Figures  21  and  22,  respectively.  The  effect  of 
a  change  in  D*  alone  over  the  range  5. 15  to  15,45  can  be  seen  from  Figures 
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14,  21,  and  22.  Comparison  of  these  figures  shows  that  the  prime  effect  of 
D*  is  to  change  both  the  dimensionless  natural  frequencies  and  the  dimen¬ 
sionless  velocity  spectral  density  at  those  natural  frequencies  present  in  the 
spectrum.  From  consideration  of  basic  vibration  theory,  this  result  is  not 
unexpected. 

It  is  again  instructive  to  examine  the  dimensionless  "peak  velocity 
spectra"  derived  from  Figures  19,  21,  and  22.  Figure  23  presents  this  com¬ 
parison.  It  is  seen  that  the  prime  effect  of  the  plate  rigidity  is  to  change  the 
lower  limits  of  applicability  of  the  "peak  velocity  spectrum"  without  signifi¬ 
cantly  changing  the  spectral  level.  This  becomes  obvious  when  it  is  recalled 
that  the  "peak  spectrum"  for  the  application  described  here  is  only  valid  above 
the  first  dimensionless  plate  natural  frequency,  and  it  can  be  shown  that  the 
first  plate  dimensionless  natural  frequency  varies  as  the  one-half  power  of  D+  . 
Comparison  of  the  lower  limiting  frequencies  in  Figure  23  bear  out  this 
relation. 

5. 1.  3  Dimensionless  Plate  Velocity  "Peak  Spectral  Density" 

As  previously  stated,  the  dimensionless  plate  velocity  "peak  spectral 
density"  is  a  useful  engineering  result  in  that  it  allows  prediction  of  the  maxi¬ 
mum  plate  velocity  spectral  levels.  Often  this  knowledge  is  sufficient  to  enable 
the  engineer  to  make  a  decision  concerning  a  design  problem  or  an  equipment 
specification.  It  would  therefore  be  desirable  to  derive  an  expression  by 
which  the  dimensionless  "peak  plate  velocity  spectrum"  could  be  predicted. 

Such  an  expression  would  have  the  further  advantage  that  it  would  show  the 
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effect  of  the  major  parameters  on  the  dimensionless  plate  velocity  spectrum 
at  the  dimensionless  plate  natural  frequencies.  The  results  of  the  above  sec¬ 
tions  provide  sufficient  data  for  such  a  derivation  over  the  range  of  param¬ 
eters  pertinent  to  the  submarine  application. 

From  equation  (4. 139)  and  the  discussion  in  Section  5. 1. 2,  it  is  obvious 


that  the  dimensionless  plate  velocity  "peak  spectral  density"  (s*,y *,<</)] 


is  proportioned  to  m25*  for  dimensionless  frequencies  above  1.932/5* . 
Section  5. 1. 2. 1  further  concludes  that,  in  the  range  1/5  <**  <  1/2,  y*  »  1/3, 
the  dimensionless  peak  spectral  density  is  independent  of  the  plate  coordinates 
*+  and  y + .  Hence, 

(5.3) 


=  (*♦> 

p  p 


1  ♦  1  ♦  1 
5*  “2  7  3 


Further,  in  Section  5. 1. 2. 2,  the  effect  of  damping  on  the  amplitude  and  slope 
of  the  dimensionless  plate  "peak  velocity  spectrum"  was  presented  in  equations 
(5. 1)  and  (5. 2).  Finally,  Section  S.  1. 2. 3  showed  that  the  effect  of  D*  was 
essentially  that  of  shifting  the  lower  dimensionless  frequency  limit  of  applica¬ 
tion  of  the  "peak  spectrum"  concept. 

By  noting  from  Figures  15,  18,  and  23,  that  the  slope  of  the  dimension¬ 
less  "peak  spectral  density"  is  constant  with  dimensionless  frequency  on  a 
semi-logarithmic  scale,  and  by  recalling  that  for  the  submarine  application 

w!,  >  L2H,  it  would  seem  reasonable  from  (5. 3)  and  the  above  discussion  to 
S* 

assume  the  following  form  for  the  plate  velocity  "peak  spectral  density": 


p 


,M  J  \  «♦“'  V,c"n> 

V‘*J 


1  ,  1  +  1 

5  =  *  =2’y  }' 

.  +  ,  1.932  . 

(i)  >  w  |  |  ,  u>  >  — ~~i —  > 


(5.4) 
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or,  in  logarithmic  form : 

10  log10  <t»;  («*)  -  Cj  ♦  20  lo|,0  M  -  20  log,0  S‘  ♦  log,0 


♦  ♦  .  1932 

"  >  wn  ’  "  >  ~ 


(5.4a) 


where  C2 ,  n ,  and  g(r*)  are  to  be  determined  from  (5. 1)  and  (5. 2). 

Using  the  values  of  M  and  8*  given  in  Section  5. 1. 2. 1,  comparison  of 
(5.1)  and  (5.4a)  yields 


-101.5-20  1o*10^~®LA=Cj+  10g/-p.Vog,o^  ♦  58.4  ♦  10  „  log,0  f  -f-  1 


from  which 


n  *=  -2 


and 


C,=  -199.9-  10  gf-L_]lo,10^. 


(5.5) 


(5.6) 


(5.7) 


Substituting  (5.6)  and  (5.7)  into  (5.4a)  results  in 
10  log104>*  (*♦)  =  -  199.9  ♦  20  log,,  M  -  20  log,0«*  -  20  log)0[-^-]+  1°  g[-f-W10 


!<,♦<! 

5  =  ”2 

.  y* 

OJ*  >  O)  j  |  , 

♦ 

CJ 

(5.4b) 


Comparison  of  (5. 4b)  and  (5.  2)  yields 


-68.5  0  log  ,0fcM  10  gi 

\ctnn/  \  cmnj 
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or 


(5.8) 


Substituting  (5.8)  into  (5.4b)  yields 

>0  1og,0«;  (a/)  « -199-9  +  20  log ,0M- 20  log, 0  «*-201og,0 


0 


1+1*1 
5*'  ■<2’y  -3* 

♦  ♦  +  1932 

01  >  tli  |  |  ,  Ol*  >  - 


«+ 


(5.9) 


or,  by  a  slight  rearrangement, 


10  log  ,„<!>♦  («/)  =  -159.9-  20  log,0 


i+i  +i 

-  <  *  <  -  ,  y-  - 

5  ”  =  2  3 
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where 


and 


«i.  ^ 


r*  2u>* 

c  m  n 

m  n 


(5.  9a) 


(5. 10) 


(5.11) 


Figure  24  compares  the  dimensionless  plate  velocity  "peak  spectrum" 
calculated  from  equation  (5.9a)  to  the  computed  spectrum  of  Figure  16. 
Agreement  between  the  dimensionless  plate  velocity  spectral  density  at  the 
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plate  natural  frequencies  and  the  dimensionless  "peak  spectrum"  are  seen  to 
be  excellent.  It  may  therefore  be  concluded  that  equation  (5.9a)  accurately 
represents  the  effect  of  major  parameters  on  the  dimensionless  plate  velocity 
spectral  density  at  the  plate  natural  frequencies  and  further  predicts  a  plate 
velocity  "peak  spectrum"  in  agreement  with  results  computed  from  equation 
(4.  139). 

5.  2  Results  of  Cavity  Acoustic  Pressure  Spectral  Density  Computations 

The  dimensionless  cavity  acoustic  pressure  spectral  density  was  com¬ 
puted  from  equation  (4. 145).  Computations  were  performed  on  the  IBM  704  and 
the  IBM  7090  digital  computers  by  iterative  processes.  Again,  the  summation 
terms  in  equation  (4. 145)  appeared  to  converge  monotonically  in  all  cases. 
However,  because  of  the  octuple  summation  in  (4. 145),  the  computational  time 
required  for  a  spectral  point  was  large,  especially  at  the  higher  plate  modes. 
Therefore,  the  solution  was  accepted  if  two  successive  iterations  agreed 
within  20  percent.  This  gave  a  computational  accuracy  to  within  +  1  decibel, 
which  compares  well  with  experimental  accuracy.  Even  with  this  fairly  large 
tolerance,  the  computation  of  the  dimensionless  cavity  acoustic  pressure  spec¬ 
tral  density  for  a  frequency  corresponding  to  the  5-1  plate  mode  required  more 
than  20  hours  on  the  IBM  704  computer.  The  IBM  7090  computer  proved  to  be 
approximately  eight  times  faster. 

Because  of  the  large  computer  time  required  per  spectrum,  the  dimen¬ 
sionless  cavity  acoustic  spectrum  was  generated  by  computing  the  spectral 
density  at  a  dimensionless  frequency  equal  to  each  dimensionless  plate 
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natural  frequency  plus  the  ratio  of  damping  to  critical  damping,  and  at  two 
points,  equally  spaced  in  frequency,  between  each  pair  of  adjacent  plate 
natural  frequencies.  It  will  be  recalled  that  the  plate  velocity  spectrum  was 
generated  similarly,  except  that  the  spectral  density  was  computed  at  three 
points,  equally  spaced  in  frequency,  between  each  pair  of  adjacent  natural 
frequencies.  This  change  results  in  a  cavity  acoustic  pressure  spectrum  less 
precise  than  the  plate  velocity  spectrum,  but  the  computer  time  requirement 
dictated  some  time  economies. 

It  will  be  recalled  that  acoustic  damping  effects  have  been  neglected  in 
the  derivation  of  equation  (4. 145),  and  that  the  only  damping  effects  are  asso¬ 
ciated  with  the  plate.  It  was  thus  felt  that  solutions  of  (4. 145)  at  or  above  the 
first  cavity  acoustic  natural  frequency  would  not  be  physically  valid.  Hence, 
computations  of  the  dimensionless  cavity  acoustic  pressure  spectrums  were 
confined,  with  two  exceptions,  below  the  first  cavity  acoustic  natural  frequency. 

To  the  author's  knowledge,  no  experimental  information  exists  for  a  case 
corresponding  to  the  theoretical  model  treated  here. 

The  purpose  of  computing  the  dimensionless  cavity  acoustic  pressure 
spectra  from  equation  (4. 145)  was  to  show  the  effects  of  major  parameters 
and  to  establish  the  range  of  validity  of  the  assumption  that  the  acoustic  back 
pressure  may  be  neglected. 

5. 2. 1  Effect  of  Major  Parameters  on  the  Cavity  Acoustic  Pressure  Spectrum 

The  data  cases  used  to  study  the  effects  of  major  parameters  were 
selected  with  the  submarine  application  in  mind.  Since  the  plate  natural 
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frequencies  and  critical  modal  damping  (b/ n  and  r *  )  are  functions  of  D* 

and  b*.  and  U*  »  1.54  throughout  this  analysis,  the  independent  major  param¬ 
eters  affecting  the  dimensionless  cavity  acoustic  pressure  spectral  density  are 

+  ♦  it  ♦  ♦  ♦  ♦  ♦ 

,  o  ,b  tr  ,D  ,c  ,d  %p  ,x  ,y  ,z  ,  uj  . 

Inspection  of  equation  (4. 145)  and  equations  (4. 131)  through  (4. 137), 

(4. 142)  and  (4. 143)  shows  that  the  effects  of  b* ,f4 , D*  ,c\d* , *  +  ,y f , z1 ,  and  a/ 
are  obscured  by  the  summation,  whereas  the  effects  of  M.S*.  and  p *  are 
explicit  in  (4. 145).  The  effect  of  the  cavity  dimensions,  b+  and  d+,  will  be 
primarily  that  of  redistributing  the  cavity  acoustic  natural  frequencies,  and 
since  these  computations  were  confined,  for  the  most  part,  below  the  first 
cavity  acoustic  natural  frequency,  the  effects  of  b+  and  d*  were  not 
investigated. 

As  previously  stated,  the  damping  (r*)  present  in  this  analysis  is  asso¬ 
ciated  with  the  plate.  Note  that  the  terms  in  equation  (4. 145)  containing  r+ 
(namely,  A  ,t  )  occur  in  exactly  the  same  combinations  and  are 

summed  over  the  same  range  of  indices  as  the  same  terms  containing  r+  in 
equation  (4. 139).  If  the  concept  of  a  "peak  spectrum"  is  again  introduced  for 
the  cavity  acoustic  pressure,  it  would  seem  reasonable  from  the  above  to  sus¬ 
pect  that  the  effect  of  damping  on  the  cavity  acoustic  "peak  spectrum"  should 


be  the  same  as  that  on  the  plate  velocity  "peak  spectrum. "  Hence,  from 
equations  (5. 1)  and  (5.  2),  it  is  assumed  that 


(5. 12) 
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and 
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(5. 13) 


where  c,  is  the  acoustic  pressure  "peak  spectral  density"  at  a  dimension¬ 
less  frequency  of  and  c2  is  the  slope  of  the  "peak  spectrum"  with 

frequency.  These  constants  will  be  evaluated  in  a  later  section. 

From  the  above,  it  remains  to  investigate  the  effects  of  Df,c*,x\y*,i*, 
and  m*. 


5.  2. 1. 1  Effect  of  Cavity  Coordinates 

Dimensionless  input  parameters  which  fall  in  the  range  of  interest  for 
submarine  sonar  applications  are 


M  =  18.9 
D*  =  10.3 

8*  =  2.31  x  10'J 

b*  =  0.6667 

r+  =  0.1 t* 

cm  n 

c*  =  454 
d*  =  0.3333 
p*  =  18.9  . 

Note  that  the  first  five  parameters  are  the  same  as  those  listed,  in  Section 
5. 1. 2. 1  for  the  plate.  This  set  of  parameters  will  hereafter  be  referred  to 
as  Case  1. 

The  effect  of  plate  coordinates  wns  investigated  for  the  following  sets  of 
cavity  coordinates  (x*,y+,z+): 


'■***H*%X. 


I 
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Because  of  computer  time  limitations,  the  frequency  range  of  interest  was 
10  <  <./  <  1000.  The  first  cavity  acoustic  dimensionless  natural  frequency  is 
1426.  The  computed  dimensionless  cavity  acoustic  pressure  spectral  densities 
for  the  Case  1  parameters  and  the  above  listed  coordinates  are  presented  in 
Figures  25  through  28,  respectively.  Close  inspection  of  these  figures  re¬ 
veals  very  small  differences  between  these  spectra  throughout  the  entire  fre¬ 
quency  range.  At  the  higher  frequencies,  it  may  be  seen  that  the  dimension¬ 
less  cavity  acoustic  pressure  spectral  density  is  slightly  greater  at  the  bottom 
of  the  cavity  (z*  =  0)  than  at  the  plate  (z*  =  -1/3).  Also  note  (from  Table  1) 
that  the  1-1,  3-1,  and  1-3  plate  modes  appear  in  the  cavity  acoustic  pressure 
spectra  of  Figures  25  through  28,  regardless  of  cavity  coordinates.  The 
above  observations  suggest  that  the  dimensionless  cavity  acoustic  pressure 
spectral  density  is  nearly  uniform  with  position  (over  the  range  of  coordi¬ 
nates  investigated)  throughout  the  cavity  at  dimensionless  frequencies  below 
i.he  fir^t  cavity  acoustic  natural  frequency.  This  is  further  borne  out  by 
Figure  29,  which  compres  the  dimensionless  cavity  acoustic  pressure  "peak 
spectra"  derived  from  Figures  25  through  28  above  the  first  plate  natural 
frequency.  It  is  seen  that  the  effect  of  cavity  position  on  the  "peak  spectra" 
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is  very  small  in  the  range  investigated;  hence,  it  will  be  assumed 

$*  (x\ y\ **....•)  «1>‘  (<->*)  i  -  x  *  -  .  y  *  -,--^*<0, 

“r  \  4  233 

.  ,  *  .  1932 

it)  '  UJ  it)  >  111 

s* 

(5. 14) 

From  the  above,  it  is  concluded  that,  in  the  range  of  space  variables 
1/3  <  x*  <  1/2,  y 4  =  1/3,  -1/3  g  z*  <  o,  the  dimensionless  cavity  acoustic  pressure 
spectral  density  is  essentially  uniform  for  dimensionless  frequencies  less  than 
2/3  of  first  cavity  acoustic  natural  frequency.  It  is  therefore  assumed  that 
the  dimensionless  cavity  acoustic  pressure  "peak  spectral  density"  is  inde¬ 
pendent  of  space  coordinates. 

5.  2. 1.  2  Effect  of  Speed  of  Sound  of  Cavity  Fluid 

A  change  of  the  cavity  fluid  involves  a  change  in  the  values  of  c*  and  p * . 
The  effect  of  p* ,  as  previously  stated,  is  explicitly  stated  in  equation  (4. 145). 
Hence,  changing  the  cavity  fluid  provides  a  simple  method  of  determining  the 
effect  of  c*  on  the  dimensionless  cavity  acoustic  pressure  spectral  density. 

To  this  end,  the  following  parameters  were  selected: 


M  =  18.9 
S’  =  0.0231 
b*  =  0.6667 
r*  =  O.lr! 

c  m  n 

D+  =  18.9 
d4  =  0.3333 

c*  =  100 
p*  0.0221  . 


This  set  of  parameters  will  hereafter  he  referred  to  as  Case  4.  A  summary 
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of  all  parameter  cases  is  presented  in  Table  2. 

TABLE  2 

PARAMETERS  USED  FOR  COMPUTATION  OF  ACOUSTIC  PRESSURE  SPECTRA 


PARAMFT  ER 

CASE  1 

CASE  2 

CASE  3 

CASE  4 

M 

18.9 

0.51 

0.0221 

18.9 

<S* 

0.0231 

0.0202 

0.0202 

0.0231 

b* 

0.6667 

0.6667 

0.6667 

0  666’ 

d* 

0.3333 

0.3333 

0.3333 

0.3333 

rVr* 

cm  n 

0.1 

0.1 

0.1 

0.1 

D* 

10.3 

0.0004 

0.0259 

10.3 

c  * 

454 

5 

5 

100 

P* 

18  9 

0.51 

0.0221 

0.0221 

wii 

1029 

0.641 

5.16 

102.9 

first  acoustic 
natural  frequency 

1426 

15  71 

15.71 

314.2 

1.932 

6* 

r 

rr\ 

00 

95.7 

95.7 

837 

The  dimensionless  cavity  acoustic  pressure  spectral  density  was  com¬ 
puted  for  the  following  dimensionless  spacial  coordinates  (*\y\z+): 
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The  frequency  range  of  interest  was  io  < «,/  <  500 .  From  Table  2,  it  is  seen 
that  this  range  includes  the  first  acoustic  natural  frequency  of  the  cavity.  The 
computed  spectra  for  Case  4  are  presented  in  Figures  30,  31,  and  32. 
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Comparison  of  Figures  30  through  32  again  shows  no  effect  due  to  cavity 
coordinates  below  a  dimensionless  frequency  equal  to  2/3  of  the  first  cavity 
acoustic  natural  frequency.  Above  this  point,  the  acoustic  pressure  spectrum 
at  the  plate  becomes  higher  than  the  spectrum  at  the  bottom  and  middle  of  the 
cavity.  This  result  is  not  in  agreement  with  the  results  of  Case  1,  which 
showed  the  cavity  acoustic  pressure  spectrum  to  be  slightly  larger  at  the  bot¬ 
tom  of  the  cavity  than  at  the  plate.  This  discrepancy  is  not  understood, 
although  it  seems  to  be  due  to  the  change  in  cavity  fluids. 

Note  further  that  the  cavity  acoustic  natural  frequency  appears  as  a  sin¬ 
gularity  in  the  spectra  of  Figures  30,  31,  and  32.  This  is  due  to  the  absence 
of  acoustic  damping. 

Comparisons  of  Figures  25  and  30,  26  and  31,  and  27  and  32  show  a 
constant  difference  in  the  dimensionless  cavity  acoustic  pressure  spectral 
density  (independent  of  frequency)  out  to  a  dimensionless  frequency  equal  to 
approximately  2/3  of  the  Case  4  cavity  acoustic  natural  frequency.  From  pre¬ 
vious  discussion  and  equation  (4. 145),  it  was  shown  that  <l>*(x\y*,z*,  <,/)  was 
proportional  to  p*\  The  above  discussion  indicates  that  the  combined  effect 
of  p*  and  c*  is  a  constant  difference  in  the  spectra  and  is  independent  of 
frequency.  It  appears  reasonable,  therefore,  to  assume 

$*<**.  y4.  z*.  <„*) 

where  n  is  an  exponent  to  be  determined.  The  difference  between  the 
dimensionless  spectra  of  Case  l  and  Case  4  is  84,9  db.  Therefore,  from  the 
above 


***** 


77 


10  io 


<T<«\  y\  «*.  .w‘)4 


H4.9 


(5.  15) 


where  the  subscripts  refer  to  the  case  number.  From  equation  (5. 15)  and  the 
values  of  Table  2, 

n  4  (5.10) 


and  thus 


<h*<i 


/) 


(5. 17) 


5.  2. 1.3  Effects  of  Damping  and  Frequency 

As  stated  in  Section  5.  2.  I,  the  only  damping  present  is  that  of  the  plate, 
and  the  damping  terms  in  equation  (4. 145)  occur  in  exactly  the  same  form  as 
in  the  equation  for  the  dimensionless  plate  velocity  spectral  density  (equation 
(4. 159)).  Therefore,  it  is  reasonable  to  assume  that  the  variation  of  the 
dimensionless  cavity  acoustic  pressure  spectral  density  with  damping  is  the 
same  as  that  of  the  dimensionless  plate  velocity  spectral  density. 

It  will  be  recalled  that  the  effect  of  damping  was  shown  by  means  of  the 
dimensionless  plate  velocity  "peak  spectra."  The  "peak  spectrum"  concept  is 
also  useful  from  an  engineering  standpoint  in  examining  the  dimensionless 
cavity  acoustic  pressure  spectra.  Because  of  the  intended  application  of  this 
work,  the  "peak  spectra"  are  only  derived  for  dimensionless  frequencies  above 
l.'Mj  iS’ .  Further,  as  the  plate  boundary  conditions  affect  the  spectral  density 
at  the  lower  modes,  and  as  the  intent  of  the  "peak  spectrum"  is  to  generalize 
the  results  into  a  useful  engineering  tool,  independent  of  plate  boundary  con- 
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ditions,  the  "peak  spectra"  only  apply  above  the  first  dimensionless  plate 
natural  frequency.  It  was  upon  the  above  arguments  that  the  assumed  forms 
for  the  damping  and  frequency  effects  on  the  dimensionless  cavity  acoustic 
pressure  '  peak  spectrum"  (equations  (!i.  12)  and  (5. 1.1))  were  based. 

From  equation  (5. 12)  and  Figure  29,  it  can  be  shown  that 

c,  -20.5  .  (5.18) 

From  equation  (5. 13)  and  Figure  29,  it  can  further  be  shown  that 

(  ,  -91.5  .  (5.  19) 


Hence,  the  effects  of  damping  and  frequency  on  the  dimensionless  cavity  acous¬ 
tic  pressure  "peak  spectrum"  at  dimensionless  frequencies  above  1.932/fi*  and 
<■>]  j  may  be  written,  taking  account  of  equation  (5. 14): 
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5.  2.  1.4  Effect  of  Plate  Rigidity 

The  plate  rigidity  (O')  originates  from  the  plate  equations  and,  like  the 
damping,  occurs  in  exactly  the  same  form  in  equation  (4. 145)  as  it  does  in 
equation  (4. 139).  Therefore,  by  the  same  arguments  used  above,  the  effect  of 
the  dimensionless  plate  rigidity  on  the  cavity  acoustic  pressure  spectral  den- 
sitv  should  be  the  same  as  its  effect  on  the  plate  velocity  spectral  density.  In 
Section  5.  1.  2.3,  it  was  shown  that  the  plate  rigidity  redistributed  the  plate 
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natural  frequencies  and  altered  the  dimensionless  plate  velocity  spectrum  at 
these  frequencies.  Its  effect  on  the  "peak  velocity  spectrum"  was  to  alter  the 
lower  limit  of  applicability  because  of  the  changes  in  <uj.  (since  the  "peak 
spectrum"  is  only  applicable  above 

The  dimensionless  cavity  acoustic  pressure  "peak  spectrum"  derived 
herein  will  have  the  same  limitations  as  the  dimensionless  plate  velocity  "peak 
spectrum."  For  the  submarine  sonar  application,  these  limitations  are  not 
severe,  but  for  aircraft  and  missile  applications,  the  cavity  acoustic  pressure 
"peak  spectra"  lower  limit  will  probably  fall  above  the  frequency  range  of 
interest;  that  is,  the  lowest  frequency  limit  of  application  of  the  "peak  spec¬ 
trum"  ( 1  .*>32  f>')  will  probably  be  near  the  upper  frequency  limit  of  interest  in 
these  applications.  In  an  effort  to  provide  some  data  in  the  range  of  interest  of 
these  applications,  and  to  better  define  the  effect  of  D*  at  frequencies  below 
1.032  5*.  dimensionless  cavity  acoustic  pressure  spectra  were  computed  for 
two  additional  data  cases.  These  were  Cases  2  and  .‘1,  and  the  parameters  used 
are  presented  in  Table  2. 

Figures  33  and  34  present  the  Case  2  results  in  the  dimensionless  fre¬ 
quency  range  o  l  <a/  s  6  for  the  dimensionless  coordinates  (1/2,  1/3,  -1/6) 
and  (1/2,  1/3,  -1/3),  respectively.  Note,  from  Table  2,  that  the  above  fre¬ 
quency  range  lies  below  1.032  6*  and  also  below  the  first  acoustic  natural 
frequency  of  the  cavity.  Again,  these  spectra  show  no  effect  of  cavity  coordi¬ 
nates  as  the  highest  frequency  falls  below  2/3  of  the  first  cavity  acoustic 
natural  frequency. 

Figures  35  and  36  present  the  Case  3  dimensionless  cavity  acoustic 
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pressure'  spectra  in  the  dimensionless  frequency  range  1.0  •  i„*  to  for  cavity 
dimensionless  coordinates  (1/2,  1/3,  -1/6)  and  (1/2,  l/.‘l,  -1/3),  respec¬ 
tively.  Again,  this  frequency  range  falls  below  I  9t2  f>‘  but  includes,  in  the 
case  of  Figure  35,  the  first  two  acoustic  natural  frequencies  of  the  cavity. 
Comparison  of  Figures  35  and  36  again  shows  no  effects  of  cavity  coordinates 
out  to  a  dimensionless  frequency  equal  to  2/3  of  the  first  cavity  acoustic 
natural  frequency. 

Table  2  shows  that  M,n‘,  and  p‘  are  the  only  parameters  varied  from 
Case  2  to  Case  3.  The  effects  of  M  and  {>'  are  known.  Hence,  Cases  2  and 
3  may  be  used  to  show  the  effect  of  D*  below  the  dimensionless  frequency 
1 .932  .  Comparison  of  Figures  33  and  34  with  Figures  35  and  36  again  shows 

the  shift  in  the  plate  natural  frequencies  with  a  change  in  D*  as  previously 
observed.  Hence,  in  order  to  show  any  frequency  dependent  effects  of  D* , 
it  would  seem  prudent  to  show  these  effects  relative  to  the  first  plate  natural 
frequency.  Knowing  the  effects  of  M  and  p' ,  and  recognizing  that  the  cavity 
acoustic  pressure  spectrum  is  independent  of  position  for  frequencies  less 
than  2/3  of  the  first  cavity  natural  frequency,  one  can  show  that,  below  this 
frequency, 


<J>*  (x\  y*,  z*.  oi*  g»*  ), 

10  log|°  . r20  lo*‘° 

«Pa  (*  ,  y  ,  z  ,  (O 


where  the  subscripts  2  and  3  refer  to  the  Case  numbers.  Figure  37  presents 
a  construction  of  f(D],  p*,  )  from  Figures  34  through  37  for  dimension¬ 

less  frequencies  below  2/3  of  the  first  acoustic  natural  frequency.  As  can  be 


seen  from  Figure  37,  the  effect  of  D*  is  a  rather  pronounced  function  of 


HI 


frequency.  A  precise  description  of  this  effect  is  difficult  from  the  above  data, 
and  the  large  amount  of  computer  time  required  to  produce  a  spectrum  pre¬ 
vented  a  more  detailed  study.  An  approximation  of  the  effect  of  !)*  at  fre¬ 
quencies  below  1  932  5*  can  be  found  by  considering  the  effect  to  be  a  constant 
equal  to  the  median  value  of  the  function  shown  in  Figure  37.  In  this  case,  one 
could  assume 


<l>  ‘  ( x  y  *,  i  * ,  (.j  ' ) 

*  ' 
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o’ 


,  I  -932  .. 

0  <  iii  <  - •  (5.  23) 
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where  m  is  an  exponent.  The  median  value  is  shown  as  a  dashed  line  in 
Figure  37,  for  which 


a  -2.48"  -2.  V 


(5.  24) 


Hence,  as  a  rough  approximation,  the  cavity  acoustic  pressure  spectral  den¬ 
sity  varies  inversely  with  the  5/2  power  of  D*  for  dimensionless  frequencies 
below  1.932/5*.  Above  this  frequency,  the  effect  of  I)'  is  most  easily  de¬ 
scribed  in  terms  of  the  dimensionless  cavity  acoustic  pressure  "peak  spec¬ 
trum,"  and  this  effect  is  to  alter  the  lower  limiting  frequency  via  altering 
In  conclusion,  a  change  in  D*  changes  all  of  the  dimensionless  natural 
frequencies  of  the  plate  and  the  associated  plate  velocity  spectral  densities. 
This  effect  is  passed  intact  to  the  cavity  acoustic  pressure  spectrum.  At 
dimensionless  frequencies  above  the  first  plate  natural  frequency  and  above 
l  932  S’ ,  a  change  in  D’  merely  changes  the  frequency  range  over  which  the 
dimensionless  "peak  spectra"  concept  applies;  that  is,  the  plate  velocity  or 
cavity  acoustic  pressure  spectral  densities  are  changed  in  such  a  way  that  the 
associated  "peak  spectra"  are  not  functions  of  o'.  However,  as  the  "peak 
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spectrum"  concept  has  been  shown  to  apply  only  above  <■  , , .  and  as  ,  is  a 
function  of  i>*.  the  lower  frequency  limit  of  application  of  the  "peak  spectrum" 
is  affected  by  l>’  At  frequencies  below  l.oy2  5' ,  the  precise  effect  of  I)' 
has  not  been  completely  investigated.  However,  it  has  been  shown  to  be  fre¬ 
quency  dependent,  as  might  be  suspected  from  the  above.  If  the  frequency 
dependency  is  neglected,  the  effect  of  n*  at  dimensionless  frequencies  below 
1.032  .S'  can  be  roughly  approximated  by 

4>|  (x*,  y',  z*,  6)*)  D*  (5.25) 

5.  2.  1.5  Dimensionless  Cavity  Acoustic  Pressure  "Peak  Spectrum" 

As  in  the  case  of  the  plate  velocity  spectrum,  the  dimensionless  cavity 
acoustic  pressure  "peak  spectral  density"  is  a  useful  engineering  result, 
allowing  prediction  of  maximum  cavity  acoustic  spectral  levels  at  dimension¬ 
less  frequencies  above  1.932/5*  and  above  the  first  plate  natural  frequency 
(01  n).  It  will  also  be  useful  to  use  the  dimensionless  cavity  acoustic  pressure 
"peak  spectrum"  in  the  next  section  to  determine  the  limits  of  applicability  of 
equation  (4. 145).  The  previous  sections  in  addition  to  equation  (4. 145)  pro¬ 
vide  sufficient  information  to  derive  an  expression  for  the  dimensionless  cavity 
acoustic  pressure  "peak  spectrum.  " 

From  Section  5.  2. 1. 1,  i  t  was  found  that  the  "peak  spectrum"  was  inde¬ 
pendent  of  cavity  coordinates  (in  the  range  of  coordinates  examined)  for  the 
dimensionless  frequency  range  between  1.932/5’  and  2/3  of  the  first  cavity 

acoustic  natural  frequency.  From  Section  5.2. 1.  2,  it  was  further  determined 

4 

that  the  "peak  spectral  density"  was  proportional  to  c  .  Further,  Section 
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5.2.  1.4  argued  that  for  frequencies  above  l  ‘M ’  <3*  the  effect  of  plate  rigidity 
was  to  shift  the  lower  limit  of  applicability  of  the  "peak  spectrum"  in  the  same 
manner  as  was  shown  for  the  plate  velocity  "peak  spectrum."  By  utilizing 
this  information  and  equation  (4. 145),  a  general  expression  for  the  dimension¬ 
less  cavity  acoustic  pressure  "peak  spectrum"  may  be  written  as  follows: 


un  =  c,/ W.v  'd’ir’.o.’) 
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(5.26) 

where  c:  is  a  constant,  and  (» * (r * . <,>’)  is  a  function  to  be  determined. 

Section  5.2. 1.3  argued  that  the  damping  (r*)  is  associated  with  the  plate 
and,  therefore,  should  enter  the  equation  for  the  acoustic  pressure  "peak 
spectrum"  in  the  same  manner  that  it  entered  the  plate  velocity  "peak  spec¬ 
trum."  Thus,  from  equations  (5.4)  and  (5.8),  it  is  assumed  that 
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In  the  above,  K  is  a  constant  to  be  determined.  The  constants  c,  and  K 
can  be  determined  from  equations  (5.20)  and  (5.21). 
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From  t’<|U.itions  (5.28)  and  (5.21), 


-•M  S  5  log,,, 


from  which 


K  -9M5  . 


(5.  29) 


(5.  50) 


Equation  (5.20)  was  hast'd  on  Case  1  variables.  Hence,  using  Case  1  param¬ 
eters,  from  equations  (5.  20),  (5.28),  and  (5.50),  one  obtains 
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from  which 


n  =  -2 


and 


10  |o„,0  C,  -250.5 


(5.51) 


(5.52) 


(5.  55) 


From  equations  (5.28),  (5.50),  (5.52),  and  (5.53),  by  a  slight  rearrangement, 
the  dimensionless  cavity  acoustic  pressure  "peak  spectral  density"  may  be 
expressed  as 


3  n  3  a  i 

(5.34) 


where  <</  is  the  first  cavity  acoustic  natural  frequency. 

*  i 

Figure  38  compares  the  dimensionless  cavity  acoustic  pressure  "peak 
spectrum"  computed  from  equation  (5.34)  to  the  dimensionless  cavity  acoustic 
pressure  spectrum  computed  from  equation  (4. 145)  for  the  Case  1  parameters 
at  the  cavity  coordinates  (1/2,  1/3,  -1/6).  The  agreement  is  very  good 
above  <./  =  150  Thus,  within  the  prescribed  limits,  equation  (5. 34)  predicts 
a  dimensionless  cavity  acoustic  pressure  "peak  spectrum"  that  agrees  well 
w.th  the  results  of  equation  (4. 145). 

5.  2.  2  Limits  of  Applicability  of  Theory 

At  the  outset  of  this  analysis,  it  was  assumed  that  the  cavity  acoustic 
pressure  was  much  less  than  the  turbulent  pressure  and  thus  provided  negli¬ 
gible  excitation  to  the  plate  compared  to  the  turbulent  pressure  excitation. 

The  calculated  spectra  and  "peak  spectra"  of  the  previous  sections  provide  the 
information  necessary  to  check  this  assumption.  Again,  it  should  be  empha¬ 
sized  that  the  calculated  cavity  acoustic  pressure  spectra" are  probably  not 
physically  valid  above  the  first  acoustic  natural  frequency. 

Horton  [43]  shows  that,  if  the  difference  between  two  random  signals  is 
10  db  (that  is,  their  ratio  is  a  factor  of  10),  the  sum  of  the  signals  is  1/2  db 
greater  than  the  larger  of  the  two  original  signals  (a  factor  of  1. 12).  Hence, 
if  the  turbulent  pressure  spectral  density  is  10  db  larger  than  the  cavity  acous¬ 
tic  pressure  spectral  density,  the  error  in  the  calculated  cavity  acoustic  pres¬ 
sure  spectral  density  would  only  be  1/2  db.  This  error  is  the  limit  of  meas¬ 
urement  accuracy.  Hence,  it  will  be  assumed  that  this  analysis  is  valid  where 
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the  cavity  acoustic  spectral  density  is  10  db  less  than  the  turbulent  pressure 
spectral  density. 

This  criterion  is  most  easily  checked  if  the  turbulent  pressure  spectral 
density  (equation  3. 1)  is  non-dimensionalized  with  respect  to  the  same  quan¬ 
tities  as  the  cavity  acoustic  pressure  spectral  density.  Let  4>’U*)  be  defined 


2.74  x  10' 5  M’S’ 

2.0  x  10’4  MJS*  J«/5 


1.932 

S’ 


1.932 
S’  ’ 


(5.  35) 


From  the  above  criterion  and  equations  (5.34)  and  (5.35),  the  limit  of 
applicability  for  dimensionless  frequencies  above  1.932/S’  and  gj,,  can 
quickly  be  established:  that  is, 
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a  comparison  of  the  Case  1  dimensionless  cavity  acoustic  pressure  "peak 
spectrum"  and  the  dimensionless  turbulent  pressure  spectrum  for  frequencies 
greater  than  1.932  /<S*.  It  can  be  seen  that  equation  (5.37)  predicts  very  closely 
the  frequency  at  which  the  cavity  acoustic  pressure  spectral  density  becomes 
10  db  lower  than  the  turbulent  pressure  spectrum. 

For  dimensionless  frequencies  below  1.932/6* ,  there  were  not  sufficient 
computations  made  to  understand  completely  the  effect  of  plate  rigidity. 
Therefore,  an  accurate  applicability  limitation  for  this  analysis  is  not  possi¬ 
ble  in  the  frequency  range.  However,  if  the  approximation  of  equation  (5.  25) 
is  used,  an  approximation  of  the  limit  of  applicability  may  be  made  for  fre¬ 
quencies  below  1.932  6*  as  follows. 

From  the  plate  velocity  spectra  computations,  regardless  of  frequency, 
it  seems  reasonable  to  assume  (as  damping  only  occurs  in  the  plate) 
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Also,  regardless  of  frequency,  from  Section  5.  2.  1.  2, 
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Further,  over  the  range  of  spacial  coordinates  investigated,  at  frequencies 
below  2  3  of  the  first  cavity  acoustic  natural  frequency, 
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By  use  of  equations  (5.38),  (5.  17),  (5.39),  and  (1.  145),  for  frequencies  below 


I .‘>32  y.  it  seems  reasonable  to  assume 
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where  c  and  f(D' ,  <■/)  are  to  be  evaluated. 
5.2.  i. 4,  f(D* ,  o')  may  be  approximated  by 

ftn'.M*)  i)'"'  2 

Therefore,  equation  (5.40)  becomes 


From  the  arguments  of  Section 
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It  remains  to  evaluate  f  , .  To  provide  a  conservative  estimate  of  the 
limit  of  application,  it  is  desirable  to  overestimate,  rather  than  underesti¬ 
mate,  the  value  of  t',.  Hence,  the  value  of  (  ,  will  be  evaluated  at  a  fre¬ 
quency  for  which  the  cavity  acoustic  pressure  spectral  density  is  a  maximum. 
From  Case  3  data  (Figures  35  and  30)  and  equation  (5.40a),  using  the  value  of 
the  spectral  density  at  one  obtains 


(  ,  -81.9 


From  Case  2  data  (Figures  33  and  34)  and  equation  (5.40a),  following  the  same 
procedure,  one  obtains 

<  ,  -80.8. 


In  order  that  the  estimate  of  the  limit  of  applicability  be  as  conservative  as 
possible,  the  largest  value  of  (  ,  was  selected.  Thus,  an  approximation  to 
the  cavity  acoustic  spectrum  for  the  frequency  range  <  >'  I  .€»32 /5*  is  given  by 
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Again,  by  use  of  the  criterion  that  the  cavity  acoustic  pressure  may  be  neg¬ 
lected  as  a  driving  force  to  the  plate  if 


,n  lo#l0 


(<,/> 


10  , 


from  equations  (5.41)  and  (5.55),  the  following  equation  must  be  satisfied  in 
order  that  equation  (4.  145)  accurately  predict  the  cavity  acoustic  pressure 
spectrum  at  dimensionless  frequencies  below  1.912  <V  : 

— L_\  .  21  log  n*  •  -26.1  .  (5.  42) 
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Using  the  applicability  limitation  of  e(|uation  (5.  42)  with  the  parameters 
for  Cases  2  and  5,  it  may  be  shown  that  Case  2  fails  the  criterion  by  approxi¬ 
mately  (il  db  whereas  Case  5  passes  by  approximately  12  db.  Figures  40  and 
41  compare  the  dimensionless  turbulent  spectrum  to  the  computed  spectrum  of 
Cases  2  and  5,  respectively,  for  dimensionless  spacial  coordinates  (1/2,  1/5, 
-1/5).  Although  the  exact  numbers  of  the  applicability  criterion  calculations 
are  not  borne  out  in  these  figures,  the  conclusions  of  the  calculations  are 
supported. 

Therefore,  it  appears  that  equation  (5.42)  may  be  used  as  an  approximate 
test  to  indicate,  for  a  particular  set  of  input  parameters,  whether  or  not 
equation  (4. 145)  will  predict  an  accurate  dimensionless  cavity  acoustic  pres- 
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sure  spectrum  at  frequencies  below  I  ')\>  o'  At  frequencies  above  l  ‘M2  .S', 
equation  (f>.  97)  should  be  used. 

It  should  also  be  noted  here,  that  the  above  criteria  must  also  be  satis¬ 
fied  for  the  calculation  of  the  plate  velocity  spectral  density  by  means  of 
equation  (4.  199)  if  the  back  of  the  plate  is  backed  by  an  acoustic  cavity.  If  the 
back  of  the  plate  is  in  vacuo  or  in  contact  with  an  unbounded  fluid  medium,  the 
above  limitations  do  not  apply. 
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SUMMARY  AND  CONCLUSIONS 

This  study  presents  an  exact  solution  for  the  acoustic  pressure  spectral 
density  and  cross  spectral  density  in  a  closed  space  behind  a  simply  supported 
plate  excited  by  boundary  layer  turbulence.  Unlike  previous  work  in  flow  in¬ 
duced  noise  in  which  models  of  the  turbulent  pressure  field  have  been  approxi¬ 
mated  by  simple  mathematical  expressions,  this  study  utilizes  an  experimen¬ 
tally  based  model  of  the  turbulent  pressure  field  as  an  input.  The  mathematical 
expression  for  this  model  is  presented  in  equation  (4.4).  The  major  assump¬ 
tion  made  in  this  analysis  was  that  the  acoustic  pressure  produced  in  the 
cavity  by  the  random,  flow  excited  vibrations  of  the  plate  could  be  neglected, 
with  respect  to  the  turbulent  pressure,  as  a  forcing  function  on  the  plate.  This 
assumption  greatly  simplified  the  mathematics  of  the  problem  but  necessitated 
an  analysis  to  determine  the  range  of  validity  of  the  results. 

Exact  mathematical  solutions  of  the  cavity  acoustic  pressure  spectral 
density  and  cross  spectral  density  were  determined  under  this  assumption, 
and  the  result^  were  presented  in  equations  (4. 126)  and  (4. 120),  respectively. 

A  useful  byproduct  of  the  mathematical  analysis  is  a  description  of  the  plate 
velocity  statistics,  which  are  presented  in  the  forms  of  the  plate  velocity 
cross  spectral  density  and  spectral  density  in  equations  (4.60)  and  (4.68), 
respectively.  It  should  be  noted  that  the  equations  for  the  plate  velocity 
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statistics  are  valid,  without  restriction,  when  the  back  of  the  plate  is  in  vacuo 
or  is  in  contact  with  a  free  acoustic  field  in  which  the  acoustic  pressure  on  the 
plate  can  be  demonstrated  to  be  negligible  with  respect  to  the  turbulent  pres¬ 
sure  excitation.  In  all  other  cases,  both  the  cavity  acoustic  pressure  and 
plate  velocity  statistics  given  in  equations  (4. 120),  (4. 126),  (4.60),  and  (4.68) 
are  subject  to  the  restrictions  imposed  by  the  above  assumption.  These  re¬ 
strictions  are  discussed  below. 

A  prime  objective  of  this  study  was  to  provide  the  submarine  sonar  de¬ 
signer  with  information  describing  the  effect  of  major  parameters  on  the 
environment  of  sonar  transducers.  The  model  used  in  this  analysis  provides 
a  fair  representation  of  the  acoustic  environment  of  submarine  sonar  trans¬ 
ducers.  In  order  that  maximum  information  be  obtained  from  selected  com¬ 
putations,  the  plate  velocity  spectral  density  and  cavity  acoustic  pressure 
spectral  density  were  rewritten  in  dimensionless  form.  The  dimensionless 
form  of  the  plate  velocity  spectral  density  is  presented  in  equation  (4. 193),  and 
the  dimensionless  cavity  acoustic  pressure  spectral  density  in  equation  (4. 145). 
Parametric  studies  and  the  range  of  validity  of  this  analysis  were  determined 
from  these  dimensionless  forms. 

Computations  of  the  plate  velocity  spectral  density,  made  with  a  digital 
computer,  compared  well  with  existing  experimental  information.  A  "peak 
spectrum,  "  constructed  from  the  computed  spectrum  by  connecting  the  major 
spectral  peaks,  proved  to  be  a  useful  concept.  From  an  engineering  stand¬ 
point,  knowledge  of  the  "peak  spectrum"  is  equivalent  to  knowledge  of  the 
maximum  expected  plate  velocity  or  cavity  acoustic  pressure  spectral  densities 
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over  some  range  of  frequencies.  This  information  is  often  sufficient  for 
making  a  design  decision  or  writing  specifications.  Comparison  of  the  dimen¬ 
sionless  plate  velocity  "peak  spectrum"  with  "peak  spectra"  constructed  from 
experimental  information  further  showed  that,  above  the  first  plate  natural 
frequency,  the  plate  velocity  "peak  spectrum"  was  independent  of  plate  bound¬ 
ary  conditions. 

From  computations  of  the  dimensionless  plate  velocity  spectral  density, 
with  input  parameters  applicable  to  submarines,  dimensionless  plate  velocity 
"peak  spectra"  were  constructed.  Analysis  of  these  "peak  spectra"  enabled 
formulation  of  a  mathematical  expression  for  the  dimensionless  plate  velocity 
"peak  spectral  density."  This  expression  is  presented  in  equation  (5.9a),  and 
has  the  advantages  that  it  shows  the  effects  of  major  parameters  on  the  plate 
velocity  spectrum  and  allows  a  quick  estimation  of  the  plate  velocity  spectrum. 
It  should  be  noted  that  the  average  time  required  for  each  plate  velocity  spec¬ 
trum,  computed  from  equation  (4. 139)  with  the  IBM  704  digital  computer,  was 
approximately  8  hours.  Thus,  the  "peak  spectrum"  concept  represents  a  sub¬ 
stantial  economic  advantage. 

The  dimensionless  cavity  acoustic  pressure  spectral  density  was  com¬ 
puted  at  various  positions  within  the  cavity  for  four  data  cases  with  a  digital 
computer.  These  four  cases,  in  addition  to  the  aforementioned  plate  data, 
allowed  determination  of  the  effects  of  major  parameters  on  the  cavity  acoustic 
pressure  spectrum  over  the  range  of  dimensionless  frequency  from  approxi¬ 
mately  i.«M2  (S'  to  2/3  of  the  first  cavity  acoustic  natural  frequency.  Above 
the  first  cavity  acoustic  frequency,  the  cavity  acoustic  pressure  spectrum 
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computed  from  equation  (4.  145)  is  not  felt  to  be  valid  since  acoustic  damping 
effects  were  neglected  in  its  derivation.  Again,  the  "peak  spectrum"  concept 
proved  useful,  and  an  equation  describing  the  dimensionless  cavity  acoustic 
"peak  spectral  density"  is  presented  in  equation  (5.  34)  over  the  range  of  di¬ 
mensionless  frequency  between  1.932/5*  to  2/3  of  the  first  acoustic  natural 
frequency. 

This  equation  describes  the  effect  of  major  parameters  on  the  cavity 
acoustic  pressure  spectrum.  The  extremely  large  amount  of  computer  time 
required  to  produce  each  cavity  acoustic  pressure  spectrum  prevented  accu¬ 
rate  description  of  the  "peak  spectrum"  over  the  entire  range  of  frequency. 

Since  the  contribution  of  the  cavity  acoustic  pressure  as  an  exciting  force 
on  the  plate  in  this  study  was  neglected,  the  results  may  only  be  used  where 
this  assumption  is  borne  out.  Comparison  of  the  turbulent  pressure  spectral 
density  and  the  cavity  acoustic  pressure  spectrum  was  made  to  determine  the 
conditions  under  which  the  results  of  this  study  are  valid.  For  dimensionless 
frequencies  above  1.932/5*.  the  condition  was  in  the  form  of  a  lower  frequency 
limit,  given  by  equation  (5.37),  below  which  the  theory  is  invalid.  For  dimen¬ 
sionless  frequencies  below  1.932/5* ,  the  condition  of  validity  could  only  be 
approximated  because  of  an  incomplete  description  of  the  effect  of  major 
parameters  in  this  frequency  range.  This  approximate  condition  is  presented 
in  equation  (5.42). 

In  summary,  the  theoretical  solution  to  this  complex  problem  and  the 
conditions  for  applicability  have  been  obtained.  From  selected  computed  plate 
velocity  and  cavity  acoustic  pressure  spectra,  dimensionless  plate  velocity 
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and  cavity  acoustic  "peak  spectra"  have  been  derived  which  show  the  effects 
of  major  parameters  on  the  plate  velocity  and  cavity  acoustic  pressure  spectra 
over  a  range  of  frequency  of  interest  to  the  submarine  sonar  designer.  The 
"peak  spectra"  have  been  shown  to  be  useful  engineering  tools  in  that  they 
allow  a  quick  estimate  of  maximum  spectral  levels  for  a  given  set  of  input 
parameters.  The  study  of  major  parameters,  specifically  the  plate  rigidity, 
at  dimensionless  frequencies  below  should  be  completed.  There  is 

also  a  need  for  extensive  experimental  work  in  this  area  to  verify  these  theo¬ 
retical  results  and  their  range  of  application. 
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Figure  2.  Longitudinal  Component  of  the  Turbulent  Wall  Pressure  Cross  Power  Spectral  Density  [33j 
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Figure  4.  The  Dependence  of  Convection  Velocity  on  Frequency  [33] 
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Figure  8.  Measured  Response  of  a  3.  5  *  3.  5  *  0.  01  inch  Steel  Plate  to  Turbulent  Boundary  Layer 

Excitation  [l5] 


Figure  9.  Computed  Response  of  a  3. 5  *  3. 5  *  0. 01  inch  Steel  Plate  to  Turbulent  Boundary  Layer 

Excitation 
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Plate  Excited  by  Boundary  Layer  Turbulence 


Figure  12.  Computed  Dimensionless  Plate  Velocity  Power  Spectrum  at  Dimensionless  Plate 
Coordinates  (1/5,  1/3);  10  Percent  Critical  Damping 


Figure  13.  Computed  Dimensionless  Plate  Velocity  Power  Spectrum  -’t  Dimensionless  Plate 
Coordinates  (1/3,  1/3);  10  Percent  Critical  Dam  g 
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Figure  15.  Effect  of  Plate  Coordinates  on  the  Dimensionless  Plate  Velocity  "Peak  Spectrum 
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Figure  16.  Computed  Dimensionless  Plate  Velocity  Power  Spectrum  at  Dimensionless  Plate 
Coordinates  (1/2,  1/3);  1  Percent  Critical  Damping 
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Figure  17.  Computed  Dimensionless  Plate  Velocity  Power  Spectrum  at  Dimensionless  Plate 
Coordinates  (1/2,  1/3);  5  Percent  Critical  Damping 
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Coordinates  (1/2,  1/3);  D*  =  5. 15 
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Figure  22.  Computed  Dimensionless  Plate  Velocity  Power  Spectrum  at  Dimensionless  Plate 

Coordinates  (1/2,  1/3);  D*  =  15.45 
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Figure  23.  Effect  of  Plate  Rigidity  on  the  Dimension! e«8  Plate  Velocity  "Peak  Spectrum 
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Figure  24.  Comparison  of  Computed  Dimensionless  Plate  Velocity  Power  Spectrum  with 

Dimensionless  Plate  Velocity  "Peak  Spectrum" 
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Figure  25.  Computed  Dimensionless  Cavity  Acoustic  Pressure  Power  Spectrum  at  Dimensionless 

Cavity  Coordinates  (1/2,  1/3,  0)  for  Case  1  Parameters 
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Figure  26.  Computed  Dimensionless  Cavity  Acoustic  Pressure  Power  Spectrum  at  Dimensionless 
Cavity  Coordinates  (1/2,  1/3,  -1/6)  for  Case  1  Parameters 
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Figure  27.  Computed  Dimensionless  Cavity  Acoustic  Pressure  Power  Spectrum  at  Dimensionless 
Cavity  Coordinates  (1/2,  1/3,  -1/3)  for  Case  1  Parameters 
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Figure  28.  Computed  Dimensionless  Cavity  Acoustic  Pressure  Power  Spectrum  at  Dimensionless 
Cavity  Coordinates  (1/3,  1/3,  -1/6)  for  Case  1  Parameters 
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Figure  35.  Computed  Dimensionless  Cavity  Acoustic  Pressure  Power  Spectrum  at  Dimensionless 
Cavity  Coordinates  (1/2,  1/3,  -1/6)  for  Case  3  Parameters 
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Figure  36.  Computed  Dimensionless  Cavity  Acoustic  Pressure  Power  Spectrum  at  Dimensionless 
Cavity  Coordinates  (1/2,  1/3,  -1/3)  for  Case  3  Parameters 


Figure  37.  Effect  of  Plate  Rigidity  on  the  Dimensionless  Cavity  Acoustic  Pressure  Power 
Spectrum  at  Dimensionless  Frequencies  Below  1.932/5* 


Figure  39.  Comparison  of  Dimensionless  Turbulent  Pressure  Power  Spectra  with  Dimensionless 
Cavity  Acoustic  Pressure  "Peak  Spectrum"  for  Case  1  Parameters 
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Figure  41.  Comparison  of  Dimensionless  Turbulent  Pressure  Power  Spectra  with  Computed 
Dimensionless  Cavity  Acoustic  Pressure  Power  Spectrum  for  Case  3  Parameters 
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II  lU^^LIMINT  ARY  NOTEI 


1}  abstract 


12  SPONSORING  MILI  T  ARY  ACTIVITY 

U.  S.  Navy 


An  analytic  solution  is  obtained  for  the  acoustic  pressure  statistics  in  a  closed 
rectangular  shaped  cavity  behind  a  simply  supported,  rectangular  plate  excited  by 
boundary  layer  turbulence.  The  contribution  of  the  cavity  acoustic  pressure  is  neg¬ 
lected  as  contributing  to  the  plate  excitation,  leaving  only  the  turbulent  pressure 
fluctuations  as  the  exciting  force.  The  mathematical  model  for  the  turbulent  pressurt 
statistics  is  based  on  that  of  Corcos,  which  agrees  well  with  experiment.  A  bypro¬ 
duct  of  this  analysis  is  an  analytic  solution  for  the  turbulent  flow  excited  plate  vibra¬ 
tion  velocity  statistics.  The  plate  velocity  and  cavity  acoustic  pressure  statistics  are 
expressed  in  the  form  of  cross  power  spectral  densities  and  power  spectral  densi¬ 
ties.  Dimensionless  forms  of  the  plate  velocity  spectral  density  and  cavity  acoustic 
pressure  spectral  density  are  developed. 


The  dimensionless  plate  velocity  spectral  density  and  dimensionless  cavity  acous¬ 
tic  pressure  spectral  density  were  computed,  by  means  of  a  digital  computer,  for 
se’ected  values  of  dimensionless  input  parameters.  From  these  computed  dimension¬ 
less  spectra,  the  effects  of  major  parameters  on  the  plate  velocity  spectral  density 
and  the  cavity  acoustic  pressure  spectral  density  were  determined. 
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